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ABSTRACT 


Investigation of the Effect of Two-Dimensional Cavities on Boundary Layers in 
an Adverse Pressure Gradient 
by 
Richard J. Margason 


The present investigation evaluated one aspect of the feasibility of the use 
of multiple cavities as an airfoil high-lift device. The effects of cavities on the 
boundary layer characteristics in several pressure gradients were determined 
experimentally and computationally. Experimentally, it was found that up to four 
cavities could be deployed with only a small change to the boundary layer 
profiles downstream of the cavities and without significantly modifying the 
resultant streamwise pressure distribution. From the computational results for 
both of the wind tunnel test section lengths used in the experimental 
investigation, it was found that a grid which provided a converged solution in 
less than a few hundred iterations was needed before a reasonable comparison 
with experimental data could obtained. It was also found for these converged 
solutions that the appropriate grid clustering and density as well as the cell size 
required for a satisfactory solution was not always apparent before comparing 
computational results with experimental data. Overall, the investigation results 
show that a multiple cavity high-lift concept may be feasible. | 
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NOMENCLATURE 


freestream speed of sound, ft/s 
test section cross-section flow area, ft 2. elemental area, ft 2 
wing aspect ratio, b’/S 

wing span, ft 

boundary condition 

experimental boundary-layer data run 
airfoil chord, ft; speed of sound, ft/s 


ee . 27 
local skin-friction coefficient, cy = =F 


pU* 


2a 


pressure coefficient, (pj-p..)/q 
local skin-friction coefficient integrated over a specified 


distance 
constant of integration in the log law equation, 


1 
U* =—In{y*)+C 
ely" }+ 
drag coefficient, D/gS 


cavity drag or friction coefficient, D/qLW 
lift coefficient, L/qS | 


computational fluid dynamics 

experimental pressure coefficient data run 

cavity depth, ft; drag force, lbf; or body diameter, ft; finite 
difference delta term a 
convection terms in the incompressible Navier-Stokes 
equations 

span efficiency, from 1/m ARe 

frequency, I/s 

fluid-in-cell method 

ratio of first universal moment to the second universal 

U 


=) 4%) 


billion floating-point-operations per second 





moment, in ( 


graphic user interface 








SJSU 





step height or tunnel height, ft — 
von Karman shape factor, 6, /9, 


incompressible two-dimensional Navier-Stokes code 
identity vector 


~ Jacobian 


speed of vortices over a cavity ratioed to freestream velocity 
ratio of convection velocity of vortices to freestream velocity 
length or mixing length, ft; local | 

length, cavity length, or run length, ft; or lift force, lbf 

mode number for the peak amplitude frequency in the cavity, — 
M,tMy, integer | 
complete wavelengths of acoustic radiation, integer 
complete wavelengths of vortex motion, integer 
Mach number, U/a 

million bytes 

million Hertz, 1/ sec 

million words 

number of cavities _ 

number of wavelengths of the fundamental frequency 
contained by the cavity length, integer 

number of grid points in a given direction — 

Strouhal number, equation 1 

pressure, Ib/ft2 

panel method - Ames Research Center 

dynamic pressure, Ibf/ft2 

random access memory 


run number; residual vector 
Reynolds number, pUx/u 


Reynolds number, pUy / wu 

autocorrelation between two velocity measurements 
separated by time T 

reference area or wing area, ft2 

tangential grid spacing scale factor 

San Jose State University 





scanivalve 
time, s | 
cycle time between arrivals of a specified phase of acoustic 
radiation at the front cavity wall lip just as a shear layer 
vortex leaves, s | | 
temperature, °F 
time-averaged local streamwise or normal velocity 
components, ft/s 
fluctuating streamwise or normal velocity components, ft/s 
friction velocity, (t/p , ft/s 
velocity, ft/s | 
phase velocity of fundamental frequency, ft/s 
non-dimensional inner layer velocity, u/u* 
Coles universal wake function, equation 32 
cavity width (spanwise in a tunnel), ft 
distance in the X, Y, or Z directions, ft 
Cartesian axis system 
normalized grid dimension in the x direction 
initial upstream y coordinate, ft 
non-dimensional inner layer distance normal to the wall, 
yu'/v 
normalized grid dimension in the y direction, y/h 
constant of proportionality in acoustic frequency relation, f = 
ft/s; angle-of-attack or diffuser angle, deg.; grid refinement 
factor, equation 36. 
Clauser equilibrium parameter, ° oe Robert's 

0 
transformation stretching parameter, equation 36; 


pseudocompressibility constant, equation 30. 
boundary-layer thickness where u = 0.99 U,, ft 


boundary-layer displacement thickness, \y (i _ fy , ft 


CO 


boundary-layer momentum thickness, ‘ os [i _ fy , ft 


XI 





uniform grid spacing before transformation to physical space 
Clauser defect thickness, i.e.: the first universal moment 
about € = 0, -3|- f(n)dy 

d phase difference between two probes, s; test section diffuser 
| angle, deg.rn grid dimension before transformation to 


a" 


physical space 

Y time lag factor to account for the time increment between the 
passage of a vortex and the emission of a sound pulse at the 
downstream cavity edge, s 


N| grid dimension before transformation to physical space 

K von Karman constant used in log-law equation, 
Ut = —In(y*)+C 

rX wavelength, ft 

Me cavity sound wave length, ft 

hy vortex spacing over a cavity, ft 

A wing sweep, deg. 

LL coefficient of viscosity, Ibf-s/ft 2 

p density, lbf-s/ft 4 

II Coles wake parameter, 0.8(B + 0.5)” 

Cn Clauser non-dimensional variables, (u-U)/u and y/d 

7 shear stress, Ibf/ft 2; time increment, s 

V kinematic viscosity, L/p, ft 2/s 

QO velocity boundary condition vector, equation 33 

W stream function 

subscripts: 

a acoustic radiation 

atm | atmospheric 

Cc | calculated 

D | diameter 

e boundary-layer edge 

g.c. cavity geometric center 

1 step 





superscript: 


n 





laminar, local, or length 


cavity length 

minimum 

maximm 

momentum 

cavity number | 
wall or flat plate; initial value or location 
pitot probe 

pressure 

reference 

Static 

tunnel surface 

shear stress 

turbulent; or total 
tunnel 

vortex motion; viscous 
cavity vortex center 
streamwise 

normal to the wall 
stations; positions 
freestream 


pseudotime step 
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l. INTRODUCTION 


Aircraft wings are usually sized by their cruise requirements which occur at 
a relatively low lift coefficient. The increased lift coefficient needed for low- 
speed flight, including take-off and landing, usually requires the deployment of 
high-lift devices. These devices increase either wing area and/or lift coefficient at 
a given angle-of-attack. Several examples of high-lift devices are shown in Figure 
1. Most modern aircraft use slats and flaps (Figure 1(a)) which extend to increase 
wing area and deflect to increase camber (i.e., lift coefficient). Additionally, less 
conventional high-lift concepts include trailing-edge tabs (Figure 1(b)), variable- 
camber airfoils (Figure 1(c)), and spanwise blowing. The pressure coefficient on 
airfoils typically has a favorable gradient from flow stagnation location near the 
leading edge to the minimum negative pressure coefficient peak which is usually 
located several percent of the chord length from the leading edge. 


For high-lift conditions, the pressure coefficient rises from the negative 
pressure peak back to the ambient static pressure near the trailing edge. This is the 
adverse pressure-gradient region (dp/dx > 0). If the pressure gradient is too large 
the lift is reduced because the flow will separate on the wing upper surface and the 
suction pressure peak is reduced and moves toward the leading edge. A cavity or 
multiple cavities could be useful in delaying separation to a more adverse pressure 
gradient. The usefulness of multiple cavities for maintaining attached flow in an 
adverse pressure gradient is determined by their effect on the viscous flow in the 
boundary layer downstream of the cavities. Maintenance of attached flow requires 
that the flow over the cavities and downstream of them retain enough momentum 
to overcome the loss due to an adverse pressure gradient, shear layer flow 
gradients, and viscous dissipation. The effects are greater near the wall because of 
the reduced velocity near the surface. If the velocity at the wall is reduced to zero 
the flow separates and produces flow forward, opposite the freestream direction. 
As a result, boundary layer assumptions are no longer valid. 


High-lift aerodynamics has been the subject of numerous investigations as 
described in a classic review by A. M. O. Smith [Ref. 1] about twenty years ago. 
This review described conventional high-lift devices and included a list of the ten 











most important basic theoretical problems of high-lift aerodynamics. The list 
emphasized development of computational methods capable. of representing three- 
‘dimensional flow with boundary-layer separation and merging multi-element 
airfoil boundary layers. There have been many experimental investigations and 
applications of computational methods to high-lift aerodynamics during the past 
twenty years. As a consequence, there has been significant progress toward the 
solution of the basic research needs listed by Smith. Additionally, the review 
indicated the need for new "inventions" such as spanwise blowing and trapped- 
vortex concepts to provide improved lifting-surface flow control. 


There have been several attempts to develop the trapped-vortex concept. 
Even Leonardo da Vinci observed and sketched very realistic recirculating eddies 
due to a bluff obstacle (or spoiler) in the fifteenth century. Recently it was 
suggested [Ref. 2] that several spoilers located at different chordwise locations 
could be deployed to form several cavities in the chordwise direction (Figure 1(d). 
This may provide a means of increasing upper surface camber and thereby 
increase wing lift. The spoilers could be deployed as a high-lift device for steady 
flow conditions. In other scenarios, the spoilers could be deployed either for 
transient maneuvers of a combat aircraft at high angles-of-attack or for increasing 
the lift of a retreating rotorcraft blade. In still another application, the sudden 
opening of a cavity near the leading edge (Figure 1(e)) of a retreating rotorcraft 
blade may suppress or delay dynamic stall. Cavity flows may also be caused by 
finned surfaces, windows, bomb bays, landing gear bays, finned heat transfer 
surfaces, and other surface imperfections. Some cavity flow research has been 
conducted specifically to study the more general problem of flow separation. 


Detailed cavity experiments were first reported in the early 1950's. Since 
then a large number of investigations of cavity flow have been conducted. This 
research has most often concentrated on a single, two-dimensional rectangular 
cavity in a uniform freestream flow with no streamwise static-pressure gradient. A 
survey of the literature was made to locate investigations of either the effect of 
cavities on boundary-layer characteristics, especially for high-lift conditions; i.e., 
in an adverse pressure- -gradient. No experimental data were found in the literature 
search which documented the effect. of either adverse pressure gradients or 
multiple cavities on the boundary-layer development or on the pressure 








distributions in or near the cavities. No data for either single or multiple cavities 
were found which documented their effect on the local flow in the vicinity of the 
cavities. These effects include the influence of a door or spoiler which closes the 
cavity opening. 


A basic understanding is needed to determine how multiple cavities could 
be used on an airfoil to either delay flow separation or to increase its camber and 
obtain an increased lift at a given angle-of-attack with suitable drag characteristics. 
The present investigation simulated a lifting airfoil using a wind tunnel with test 
section liners to represent different pressure-gradients which could be associated 
with an airfoil at selected angles-of-attack. This investigation concentrated on the 
effect of cavity flow on the attached surface boundary-layer characteristics in 
adverse pressure-gradients. Conditions where flow separation was present were 
not considered. The purpose of the present investigation was to experimentally 
and computationally evaluate the effect of 1, 2, or 4 two-dimensional spanwise 
cavities with square cross-sections on the boundary-layer characteristics in both a 
constant pressure flow and for two adverse pressure gradient flows. The 
experimental investigation was conducted in the San Jose State University (SJSU) 
12" by 12" low speed tunnel. Computational results were also obtained using a 
numerical solution to the incompressible Navier-Stokes equations. 











(a) Flap and slat with chord extension. 


(b) Tabs on airfoil and flap trailing edges. 





(c) Variable camber airfoil. 


(d) Multiple spoilers used to form multiple cavity concept. 
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_(e) Leading edge cavities to delay dynamic stall. 


Figure 1. High lift concepts. 











If. PREVIOUS CAVITY FLOW RESEARCH 


Two-dimensional boundary layers in subsonic flow are well defined as 
described in a very complete survey paper by Ligrani [Ref. 3] and textbooks such 
as Schlichting [Ref. 4], White [Ref. 5], and Cebeci and Smith [Ref. 6]. At low 
Reynolds numbers, i.e. short run lengths, in low-turbulence flow the boundary 
layer is usually laminar. For high-lift applications, the boundary layer rapidly 
transitions to a fully turbulent boundary layer whose mean velocity profile may be 
represented by the 1/7 power law (u/Uco = (y/d)!//). Additional work on 
boundary layers and the effects of skin friction, adverse pressure gradient, surface 
roughness, and shear-layer flow are described in References 17 through 31. 


Nearly 40 years ago, Roshko [Ref. 32] and Krishnamurty [Ref. 33] 
conducted two of the first relatively-detailed investigations of cavity flow. These 
experiments did a good job of identifying the important flow characteristics 
associated with a single cavity. In the present paper, the cavity streamwise 
dimension, length, is designated as L; the dimension normal to the freestream, 
depth, is D; and the distance across the tunnel and normal to the flow direction, 
width, is W. In general, cavity flows exhibit various steady and unsteady 
phenomena. The upstream boundary layer separates at the cavity lip to form a 
shear layer over the cavity. The shear layer then reattaches (1) either on the cavity 
floor or (2) on the downstream cavity wall or downstream of the cavity. The 
cavity is considered to be closed if the shear layer reattaches on the cavity floor 
and then recirculates within the upstream end of the cavity. Typically this flow 
can occur when L/D > 4. The cavity is considered to be open when the shear layer 
reattaches near the rear lip or downstream of the cavity. This flow is typical of 
deeper cavities where L/D < 1. For L/D values between 1 and 4, the d5/L strongly 
influences whether a cavity is open or closed. At certain conditions disturbance 
waves are fed back to the upstream cavity lip, the original disturbance source. 
This feedback loop can amplify the disturbance waves and create large oscillating 
pressure waves and noise. In these circumstances the unsteady flow can dominate 
the cavity flowfield. 


A. MEAN FLOW CHARACTERISTICS 


Roshko [Ref. 32] used a cavity length L of 4 inches and systematically 
varied the cavity depth D to study the time-averaged flow at low Mach numbers 
(0.06 to 0.20) for D/L ranging from 0.02 to 2.50. Pressure coefficients on the 
cavity walls and floor were measured and the skin-friction coefficients were 
calculated. The structure of the flow in the cavity was also observed and found to 
be a function of the cavity depth-to-length ratio D/L. For very shallow cavities 
(D/L < 0.1) the shear layer above the cavity reattached to the cavity floor forming 
a recirculation region on the upstream corner of the cavity. As the varying cavity 
D/L approached 1, the shear layer attachment location moved aft along the cavity 
floor to the rear-wall intersection and up the wall to near its intersection with the 
freestream surface. For a cavity with a square cross section (D/L ~ 1) there was 
steady flow in a single predominant vortex with secondary vortices in the corners 
of the cavity. For D/L from 1 to 2.5, a single vortex continued to dominate the 
cavity flow. 


It was concluded that the drag due to the cavity was almost entirely due to 
the pressure on the cavity walls. The drag contribution from changes in the skin 
friction coefficient on the cavity surfaces was relatively small. Further, the friction 
forces on the cavity walls were found to be small compared with the pressure 
forces. The drag is analyzed in detail in ssection II.C. While the friction forces 
due to the cavity vortex were negligible compared to the pressure drag, they do 
play a role in determining the vortex stability. An indication of this stability is 
shown in Figure 2(a) by the variation of the pressure coefficient measured at the 
middle of the cavity floor as the cavity depth was increased from D/L of 0.02 to 
2.50. Steady pressures were measured when D/L was less than 0.50 and when D/L 
was between 0.87 and 2.00. For D/L from 0.50 to 0.87 and D/L > 
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(a) Variation of pressure coefficient measured at the middle of the cavity 
floor [Ref. 32]. 
Figure 2. Effect of cavity D/L ratio on local pressure-coefficients. 








2.00 the pressures were unsteady. The change from unstable flow to stable flow as 
D/L was varied through 0.87 was quite distinct. There were two stable states for 
this range of D/L. When the cavity vortex was steady (D/L < 0.50 and 0.87 < D/L 
< 2.00), the flow over the cavity was also steady. 


A single, stable vortex was formed by the deflection of a portion of the 
shear layer at the downstream cavity edge into the square cross-section cavity. 
The relatively high pressure on the cavity wall in that vicinity accounts for most of 
the drag. The pitot pressure at the top back corner is shown in Figure 2(b). This 
pressure measurement is analogous to a Preston-tube measurement and it is related 
to the local skin-friction coefficient. It is intended to give a measure of the 
_ pressure near the top of the rear cavity wall. There is an unsteadiness shown for 
D/L between 0.50 and 0.87. Above D/L of 1 there is a hysteresis which shows that 
the pressure coefficient is dependent upon whether D is increasing or decreasing. 
For the square cavity (L/D = 1), the skin-friction coefficient was calculated from 
the boundary-layer profiles with the following results: 


Boundary-Layer Profile Location 
I - 0.375 L upstream of the cavity front wall 0.0015 
If - 0.125 L downstream of the cavity rear wall 0.0011 


Ill - 0.375 L downstream of the cavity rear wall 0.0012 


IV - 0.375 L downstream of the cavity rear wall 
location without the cavity 1 










Another investigation by Maull and East [Ref. 34] studied cavities at 
low subsonic velocities using oil-flow and surface static-pressure distributions. 
They found that the flow steadiness depended on cavity width as well as 
streamwise length and depth. For a cavity width-to-length ratio W/L of 9, a non- 
uniform spanwise variation of pressure coefficient was most notable below D/L of 





0.85 and between D/L of 1.5 to over 2. These regions are roughly consistent with 
the steadiness of Roshko's pressure coefficient data in Figure 2. There were 
regions of uniform pressure distribution near D/L of 1 and 2.5. There was a very 
rapid switch from non-uniform pressure-coefficient variation for D/L = 0.85 to a 
nearly uniform pressure-coefficient distribution for D/L = 0.86. These results 


appeared to also depend on the ratio of boundary-layer thickness to cavity length 
O/L. 


Rossiter [Ref. 35] investigated the flow over rectangular cavities at 
subsonic and transonic speeds. The investigation was conducted in a 2' by 1.5 
transonic tunnel. The cavity had a 4" width, an 8" length, and a depth which was 
varied from 0.8 to 8.0 inches. Steady and unsteady pressures were recorded on the 
cavity floor and downstream of the cavity. While the flow was highly unsteady, it 
was useful to briefly discuss the nature of the time-average flow over cavities. For 
very shallow cavities, the flow over the front and rear walls may be considered 
independently as the flow down and up a step, respectively. The airflow will 
separate from the front edge and reattach at some point along the floor of the 
cavity. The pressure in the separated region will be lower than the freestream 
pressure. This occurs because the freestream flow speeds up as it enters the cavity. 
Then the pressure rises at the attachment point. As the airflow approaches the rear 
wall, it is slowed and then the pressure increases until a position is reached where 
the boundary layer separates to form the boundary layer ahead of the rear wall. 
The boundary layer will usually reattach at a location downstream of the cavity. 
As the depth-to-length ratio of the cavity increases, the attachment and separation 
points on its floor will move closer together until a reverse flow develops between 
the high pressure region ahead of the rear wall and the low pressure region behind 
the front wall. A large vortex then forms within the cavity. 


As shown by Roshko [Ref. 32] the mean flow pattern depends on the . 
length-to-depth ratio L/D of the cavity. Rossiter compared his data with Roshko's 
as shown in Figure 3. The results are not directly comparable between these tests 
because of different length-to-depth ratios and the relatively thicker boundary 
_ layer approaching the cavity in Rossiter's investigation. For the very shallow 
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(b) Pressure coefficient measured at the top back corner of a cavity for a 
range of depth-to-length ratios [Ref. 32]. 
Figure 2. Concluded. 
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cavities (L/D > 8) the two pressure rises associated with flow attachment and with 
separation on the floor of the cavities may be seen. Immediately behind the cavity 
front wall the pressure reaches a low value. The pressure increases from the center 
of the closed vortex in the corner at the front wall to the rear end of the vortex. 
The pressure levels off and then increases again as the downstream wall is. 
approached. At a length-to-depth ratio of 8, the two pressure rises have merged 
and extend from x/L ~ 0.4 to 0.9. Ata length-to-depth ratio of 6, the pressure 1s 


Rossiter [Ref. 44] 
— —-— Roshko [Ref. 41] 





Figure 3. Comparison of Pressure coefficients measured on cavity floors and on 
the surface downstream of the cavity by Rossiter [Ref. 35] and Roshko [Ref. 32]. 


almost constant (Cp ~ -0.02) along the floor of the cavity indicating that the two 


separations have combined. For the deeper cavities (L/D of 4 and 2) there is a 
decrease in the pressure on the floor associated with the high airspeeds at the 
periphery of the cavity vortex. As Mach number is increased, the flow attachment 
point on the floor of the shallower cavities moves downstream so that the cavities 
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become effectively deeper. As a result, the pressure distributions for the shallower 
cavities show a large variation with Mach number whereas the pressure 
distributions in the deeper cavities are comparatively independent of Mach 
number. 


Other investigations, such as Charwat et al [Ref. 36 and 37], have studied 
Cavities in supersonic flow and found that the same general flows exist as those 
found at subsonic speeds. Three-dimensional cavities which relate to landing gear 
or weapons bays were studied by Plentovich [Ref. 38]. Greater unsteadiness was 
found at the highest Reynolds number and the pressure distributions were found to 
be sensitive to the thickness of the boundary layer entering the cavity. Another 
application of cavities is found in aircraft where telescopes are used for various 
observations. The objective of the investigation by Buell [Ref. 39] was to 
minimize the shear-flow disturbance using an antiresonance device. Devices were 
developed which suppressed cavity shear-flow resonance and produced thinner 
shear layers. This shear-layer flow suppression enabled better light transmission 
to the telescope. However, the suppression of resonance also made the boundary 
layer downstream of the cavity more susceptible to flow separation. | 


B. FLOW INDUCED OSCILLATIONS 


There have been reviews by Rockwell and Naudascher [Ref. 40] and by © 
-Komerath, Ahufa, and Chambers [Ref. 41] which emphasize cavity flow-induced 
oscillations. Both papers used the classification of self-sustaining flow oscillation 
over cavities first presented by Rockwell and Naudascher and reproduced in 
Figure 4 from reference 40. Three flow-interaction categories were identified: (1) 
fluid dynamic, (2) fluid resonant, and (3) fluid elastic. In many situations more 
than one of these interactions may be involved. The fluid-dynamic interactions 
involve coupling between oscillations of the shear layer over the cavity with the 
flow inside the cavity. Many of the oscillating flows at low speeds over shallow 
cavities (L/D > 1) fall under this category. The mechanisms involved are believed 
to arise from shear-layer instability and vortex shedding. Large-scale coherent 
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structures present in the shear layer are known to play the major role in such 
interactions. Fluid-resonant interactions are flow oscillations which appear to be 
controlled by the acoustic modes of the cavity. These are usually encountered in 
cavities which have large depths normal to the flow direction (L/D < 1). These 
fluid-resonant oscillations are observed in flows over cavities at high Mach 
numbers. Fluid-elastic flows involve interactions between the shear layer over the 
cavity and the elastic boundaries of the cavity. These interactions may cause 
sonic-fatique problems. 
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Figure 4. Classification of cavity flows by Rockwell and Naudascher [Ref. 40]. 


Krishnamurty [Ref. 33] varied the rectangular cavity L/D ratio in a study of 
the sound radiated by cavity flow. The cavity had a constant depth of 0.1" and 
was located in a flat plate. A trip wire and the flat-plate angle-of-attack were 
varied to change the boundary layer from laminar to turbulent. The cavity length 
was varied from 0 to 2 inches. The investigation was conducted in a blowdown 


13 








wind tunnel over a Mach number range from 0.4 to 1.5. Data included hot wire 
measurements of the boundary layer, magnitude and directionality of the radiated 
sound, as well as, schlieren and interferometer flow visualization. Below a - 
minimum cavity length it was found that the shear layer flowed over the cavity and 
reattached downstream of the rear cavity wall. For these conditions there were no 
acoustic oscillations. 


Above this cavity length, the acoustic frequency was inversely proportional 
to the cavity length f= a/L. The constant of proportionality was different 
for laminar and turbulent boundary layers. This dependence is presented in Figure 
5 for both laminar and turbulent boundary layers. While there was a dominant 
frequency for laminar flow, there were two frequencies, low and high, for the 
turbulent case. Krishnamurty used Strouhal number to obtain a dimensionless 
frequency. If the freestream velocity is used as the characteristic velocity, the 
Strouhal number becomes 


No = fL/U,, =a/U,, | (1) 


At low Mach numbers there was little sound radiation directionality; at increased 
Mach numbers the sound became more intense and directional. The radiation 
pattern was observed using a schlieren system. It was shown that the unsteady 
pressure oscillations were greater in a laminar boundary layer than in a turbulent 
boundary layer. 


Rossiter [Ref. 35] found that the unsteady pressures contain both random 
and periodic components. The random component predominates in the shallower 
cavities (L/D > 4) and was most intense near the rear wall. The unsteady pressure 
had a smooth spectrum over a broad frequency band showing the random character 
in shallow cavities. For very shallow cavities, a local region of intense pressure 
fluctuations was found where the shear layer flow attaches to the floor of the 
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(a) Proportionality constant a. 


Figure 5. Effect of Mach number on the acoustic field due to a cavity. 
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(b) Strouhal number Nog. 
Figure 5. Concluded. 


15 





cavity. There were usually two peaks of almost equal magnitude. The periodic 
component predominates in the deeper cavities (L/D < 4) and may form standing 
wave patterns with one peak whose magnitude is much larger than any other. It is 
suggested that the periodic component is due to an acoustic resonance within the 
cavity excited by a phenomenon similar to that causing edge-tones. The periodic 
pressure fluctuations may be very large. Mean values up to 0.35 times the 
freestream dynamic-pressure were measured. These results indicated the 
predominance of the periodic component over the random oscillation. Increasing 
Mach number caused an increased periodic component. The periodic component 
was reduced in the thicker boundary-layer. This suggests that the lack of periodic 
pressure fluctuations in shallow cavities is partially due to a large ratio of 
boundary layer thickness to cavity depth 6/D. This result suggested a means to 
reduce the magnitude of the pressure fluctuations. The suppression was achieved 
using a small spoiler located upstream of the cavity to effectively increase the 
boundary-layer thickness. | | 


In this investigation, Rossiter identified experimentally the principal 
characteristics of the periodic pressure fluctuations as the following: 


(1) The pressure fluctuations may contain a number of periodic 
components or modes which occur at specific frequencies and are 
designated as m = 1, 2, 3,... 


(2) The frequency of any component is inversely proportional to the cavity 
length and increases with freestream velocity. Dimensionally this suggests 


U 
(=) = g(M,R,) (2) 
L 
Over the Reynolds number range investigated it was found that the effect of 
cavity length was small. 


(3) Shadowgraphs showed that the pressure fluctuations are accompanied 
by the periodic shedding of vortices from the front lip of the cavity while 
the principal acoustic source is close to the rear lip of the cavity. 
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(4) Under some conditions, standing wave patterns occur in the cavity 
which are probably due to an acoustic resonance within the cavity. 


The regular shedding of vortices found in flow over cavities is a feature of 
the edge tone phenomenon as well as the strong acoustic-radiation. A connection 
was assumed between these two acoustic features. Further it was assumed that the 
acoustic radiation initiates the vortex shedding and that the passage of the vortices 
over the rear lip of the cavity is responsible for the acoustic radiation. If the 
average speed of the vortices over the cavity is K times the freestream velocity and 
sound waves travel upstream in the cavity at a mean velocity c then 
KU. c 


aaa (3) 





f= 
Vv 
where A,, was assumed to equal A. The phase relation between the two motions is 
unknown. 


A solution was found by identifying the particular phase of the acoustic 


radiation which reaches the front lip of the cavity at the instant that a vortex is 
shed. It is assumed that a vortex is y A, behind the rear lip-when this particular 


_ phase of the acoustic radiation leaves the source at the rear lip. It was assumed 
_ that there were m,, complete wavelengths of the vortex motion and m, complete 


wavelengths of acoustic radiation. Further it was assumed that there is a time 
interval t' between arrivals of an identified phase of the acoustic radiation at the 


front lip just as a vortex is shed. Therefore the vortex pattern moves downstream a 
distance KU,,t' in this time interval so that 


m,A, =L+yA,+KU,t (4) 
Also in time t' the internal wave system has moved a distance ct’ so that 
L=m,A, +ct' (5) 


Then t' is eliminated between equations 4 and 5. The freestream velocity U,, was 


replaced by Ma; further, it was assumed the acoustic velocities internal, c, and 
external, a, to the cavity were equal; and that the total number of wavelengths, My 


+m, , equal the mode number, m. 
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These assumptions led to a physical model of the flow which may be 
represented by the empirical frequency equation 


a 
K 
where K equals 0.57, m is an integer, and equals 0.25. When the frequency of 
_ one of the components (m-y) is close to the natural frequency of the volume of air 


in the cavity, resonance occurs. Both cavity length and depth influence which 
mode number is dominant for a particular cavity geometry. 


Additional investigations of flow induced cavity pressure 
oscillations and acoustic resonance to refine Rossiter's work were conducted by 
East [Ref. 42], Heller et al [Ref. 43], Bilanin and Covert [Ref. 44], Miles [Ref. 45], 
Block [Ref.46], Sarohia [Ref. 47], Yu [Ref.48], and Gharib and Roshko [Ref. 49]. 
The investigation by Tam and Block [Ref. 50] obtained experimental data for the 
tone-frequency characteristics for flow Mach numbers less than 0.4. It was 
observed that there was a transition from the normal mode resonance mechanism > 
to a feedback instability mechanism for disctete tone generation as flow Mach 
number increases. A mathematical model of acoustic-feedback oscillations was 
developed which, in contrast to the Bilanin and Covert model, accounted for the 
shear-layer thickness. Compared with the Rossiter model, the Tam and Block 


model added two parameters; (1) the cavity length-to-depth ratio L/D; and (2) the 
ratio of the momentum thickness of the shear layer to the cavity length 65/L. 


The acoustic-wave generation process sketched in Figure 6 (from Ref. 51) 
for a supersonic freestream was adopted. Flow-induced cavity oscillations are 
caused by the interaction of the free shear layer and the complex internal cavity 
wavetrains. The shear layer oscillated up and down near the trailing edge of the 
cavity. During the upward motion of the cycle, the fluid of the shear layer shields 
the trailing edge of the cavity from the external flow and the predominant flow is 
over the cavity with no pressure waves. When the shear layer deflects downward , 
there is flow of the freestream into the cavity which causes a transient high 
pressure region near the cavity trailing edge, which forces propagation of a 
compression wave in all directions. The convection effect of the freesteam 
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modifies the shape of the wave front as it radiates away from the cavity trailing- 
edge. | 


In Figure 6, the essential features of a typical oscillation cycle are divided 
into six time-sequential parts, identified by the letters A through F. The external- 
flow Mach waves represent a freestream Mach number of 1.5. Each part sketched 
in Figure 6 should be viewed as a typical phase of the oscillation cycle. The 
direction of the motion of each wave is indicated by an arrow. The exact wave 
representation depends on the cavity geometry (L, D, and W), the external Mach 
number, the boundary layer characteristics at the upstream end of the cavity, and. 
the freestream turbulence. The starting point for the cycle is arbitrarily selected. It 
is helpful to review the entire cycle to gain a good understanding of cavity flow 
phenomena. This cycle is now discussed starting at A: 


(A) The pressure wave from the previous trailing-edge disturbance reaches 
the front of the cavity and reflects. Another wave, moving from the front wall 
approaches the rear wall. The shear layer is above the cavity trailing edge, so the 
external flow does not produce a disturbance at the cavity trailing edge. Some 
fluid leaves the cavity at the rear. 


(B) The shear-layer waveform travels rearward and reduces the height of 
the shear layer above the trailing edge. A new compression wave begins to form 
at the rear as the flow interacts with the trailing edge and fluid is added to the 
cavity. The front compression wave has reflected off the front wall and moves 
rearward nearly in phase with the shear layer displacement. The previous 
rearward wave has reached the trailing edge. 


(C) The wave reflected off the front wall continues to move rearward in 
phase with the shear-layer displacement. The shear layer, which is now below the 
trailing edge at the rear of the cavity, forms a new forward-traveling compression 
wave as the external flow impinges on the back of the cavity. 


(D) The newly-generated, forward-traveling compression wave and the 
reflected, rearward-traveling compression wave meet and interact near the cavity 
center. 
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Figure 6. Typical oscillation cycle for cavity flow [from Ref. 51]. 
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(E) After interaction, the waves continue in their respective directions. The 
external part of the forward-traveling wave moves into the supersonic flow, thus 
causing it to be tipped more than the external flow angle. The rearward wave 
moves in the same direction as the external flow and travels at a subsonic speed 
relative to it. This subsonic relative speed explains why the rearward-traveling 
wave stops at the shear layer. At the rear, the shear layer reaches the trailing edge 
height. 


(F) The shear layer is now above the trailing-edge height. The wave 
generated at the trailing edge approaches the front of the cavity, and the reflected 
wave nears the rear of the cavity. The next step is the same as (A), and oe 
oscillation cycle repeats. 7 


For subsonic flow, the process is essentially the same, particularly as regards the 
internal wave structure. The forward-traveling wave will still be supersonic with 
respect to the external flow. The external wave structure will usually be non- 
existent. Instead, the shear layer will tend to roll into transverse vortices with the 
number dependent primarily on the cavity geometry (L/D)) and on the freestream 
Mach number. Typically there will be either two (mode II) or three (mode 
vortical structures present. 


For the purpose of computing the phases of the acoustic waves generated at 
the trailing edge of the cavity, Tam and Block [Ref. 50] made the size of the noise 
source very small. A periodic line source was used to produce this flow pattern. 
Inside the cavity the mean velocity was zero. The mathematical model of the 
effect of the interaction of the acoustic-wave field on the instabilities of the shear _ 
layer was assumed to be convective. A mean shear-layer momentum thickness 
was used in the model. The resultant equations, which are fairly complex, are 
presented in Reference 50. This model provided good agreement between 
predicted discrete-tone frequencies and Rossiter's data (0.4 < M < 1.2) and Tam & 
Block's data (M < 0.2). For the very low Mach numbers, M < 0.2, it appeared that. 
the tones were generated by the normal-mode resonance mechanism. The 
transition between the feedback mechanism and the normal-mode resonance was 
rather gradual. This suggested that the Tam and Block model may provide the 
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basis for a unified model of the flow-induced tone phenomenon to include cai 
tone frequencies throughout the Mach number nee 


Ziada and Rockwell [Ref. 52] found that impingement of mixing layers on 
solid boundaries enhanced the organization of all harmonic modes in a shear layer. 
This effect was investigated by Rockwell and Knisely [Ref. 53] specifically for 
cavity flow using a laser velocimeter. As shown in Figure 7, inserting a 
downstream cavity-impingement edge makes a large change in the unsteady nature 
of the shear flow. Autocorrelations were taken along the edge of the shear layer 
where u/U,, = 0.95 using an average of six sequential samples. These results are 


presented in Figure 7; note that the ordinate scale changes for each curve. Both 
data sets were taken at the same shear layer location of x/d, = 134. For the 
backward facing cavity, there were large variations in the autocorrelations. The 
typical period of oscillation is quite irregular and tends to increase with | 
downstream location. For the cavity flow, the enhanced organization of the flow 
is quite dramatic. In addition, there are no apparent variations in the oscillation 
frequency. This enhancement of the shear layer organization extended along the 
entire length of the shear layer. This finding substantiates the concept of 
disturbance feedback. The perturbations propagated upstream from the 
impingement surface affected the amplitude of the locally-induced pressure 
fluctuations and the resultant shear layer deformation. 


A few examples of the many additional cavity-flow investigations are cited 
below to illustrate some results pertinent to the present investigation. Franke and 
Carr [Ref. 54] investigated the effect of geometry on modification of flow-induced 
oscillations for open cavity flow in the Mach number range from 0.6 to 3.3. It was 
found that ramps at the leading and trailing edges of the rectangular cavities could 
be effective in reducing pressure amplitudes under some conditions. Some of the 
cavity configurations studied by Franke and Carr included two cavities in tandem. 
Tandem cavities were also investigated by Zhang and Edwards [Ref. 55] in 
supersonic flow. They examined both the mean and unsteady flowfields. It was 
found that two cavities whose L/D was either 1 or 3 had little effect on each other. 
In contrast, a L/D = 3 cavity completely altered the flow in a downstream L/D = | 
cavity. 
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Betts [Ref. 56] experimentally evaluated the effect of slotted walls on flow 
ina rectangular water-channel. It was found that the effect of geometry on flow 
oscillations could be modeled by empirical relations. Sarno and Franke [Ref. 57] 
evaluated several devices located upstream of a cavity which attempted to reduce 
flow oscillation. Like Buell [Ref. 39] it was found that stationary fences at the 
cavity leading edge were the most effecive flow-suppression device since they 
effectively reduced L/do. 
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Figure 7. Comparison of autocorrelations taken with and without the cavity 
impingement edge at corresponding locations in the cavity shear-layer [Ref. 62}. 








C. CAVITY DRAG 


The drag of cavities, holes, and gaps were discussed by Hoerner [Ref. 13] 


in his book on fluid-dynamic drag. Cavity drag is relatively small compared to the 
overall aircraft drag. The cavity drag coefficient Cp was typically less than 0.01 


based on a reference area equal to the cavity opening, LW. For comparison, 
complete aircraft drag coefficient is usually based on wing area; the cavity opening 


area LW of landing gear or a stores cavity is usually at least one or two orders of 
magnitude less than wing area so that the Cp contribution to aircraft drag 


coefficient is of the order of 0.001 to 0.0001 or less. The shape of the edges of the 


cavities have a major influence on the magnitude of the drag as shown in Figure 8. 
For square cavities, L/D = 1 such as those used in the present investigation, Cr 


was 0.0083. When the slope on the top of the rear wall was changed to -5°, Cp 
was reduced to 0.0072. Rounding the rear slope reduced Cp to 0.0060. The 


other sketches in Figure 8 show a cavity with the same opening area but a larger 
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Figure 8. Cavity drag coefficient Cp measured for several edge shapes [Ref. 13]. 
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cavity volume with variations in the slopes into and out of the cavity. It is shown 
that the value of Cp can be varied from 0.0070 to as high as 0.0400 by changing 
these slopes. Cavities with front slopes which direct freestream flow into the 
cavity and with rear wall slopes which inhibit flow out of the cavity greatly 
increase the drag. For minimum drag, the upstream edge should be nearly parallel 
with the freestream and the downstream edge should be bent at a small angle into 
the cavity. 


The maximum shear stress increases with the cavity L/D. Gharib and 
Roshko [Ref. 49] measured the maximum cg as 0.013. This cg value compares 
closely with 0.0125 measured by Liepmann and Laufer [Ref. 58] for a two- 


dimensional shear layer and 0.012 measured by Kistler and Tan [Ref. 59] for a 
two-dimensional cavity shear layer. The maxima of cavity c¢ profiles did not 


have a constant value. In contrast, ce was found to be constant in two- 


dimensional, self-similar free turbulent shear layers. 


Cavity drag is defined as the net force in the flow direction experienced by 
the cavity. The cavity control volume which can be used for the momentum 


balance is shown in Figure 9(a). Surface 1 is the open interface between the cavity 
and freestream flows where drag equals the momentum flux drag Dy occurs. On 





(a) Schematic of the cavity and its control volume. 


Figure 9. Cavity drag and effect of cavity length on maximum Cp. 
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On the upstream and downstream cavity walls (surfaces 4 and 2, respectively) the 
pressures are integrated to give a pressure drag Dp. On the cavity floor (surface 


3) the shear stress gives a third drag force contribution Dg. These terms are 
summed to zero for the control volume. The drag force on walls 2, 3, and 4 is 
given by | 


D=Dp+Dsg =Dyy (7) 


The cavity drag is expressed as either the direct forces on the cavity walls and 
floor or the cavity momentum flux. The drag terms in equation 7 may be 
expanded to provide the cavity momentum-drag balance in integral form as 


ee | 
D= j PA - j Pda + hy T,dA3 = f,[ -P0¥- pu' v a ha (8) 


Roshko [Ref. 32] found that the dynamic pressure in the recirculating cavity 
flow is low. On their axisymmetric experimental cavity model, which had no 
cavity floor at the centerline, there was a standing axisymmetric vortex. At the 
body centerline the vortex velocity components tended to cancel and the resultant 
velocity was nearly zero. As a consequence, the resultant drag due to shear stress, 
Dg , was negligible compared to the pressure drag, Dp. This simplified the drag 


equation to 
D = Dp=Dy 9) 


This drag force can be nondimensionalized as a drag coefficient using the 
freestream dynamic pressure and an arbitrary reference area so that 


D 


Cp =——5— = Cp, =Cp (10) 
ly PUL" Steg : 2 


Alternatively, the reference area may be the cavity opening S,e¢ = LW to obtain 
an equivalent average friction coefficient due to the cavity (or cavity drag 


D 


Cr == = Ce =C 11 
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coefficient). Both coefficients, Cp and Cp, give the same drag value. Gharib and 
Roshko [Ref. 59] choose to associate Cp with the forces on the cavity walls and 
Cpr with cavity momentum flux. For cavities in the non-oscillation mode Cp was 
near 0.0001; in the self-sustained mode, the Cr was less than 0.01; and while in 
the wake mode (L/d9 > 155 or L/D > 1.25) the Cp was between 0.06 and 0.08. 
For comparison, the estimated Cp based on boundary layer friction in the absence 


of a cavity was about 0.005. 


Measurement of the pressure distribution on the solid surfaces of the cavity 
provides an understanding of drag generation. In short cavities with non- 
oscillating flow, the pressure coefficient Cp on both the upstream and downstream 
walls tends to be slightly negative and uniform indicating essentially no drag 
contribution. In cavities with flow oscillations (84 < L/39 < 155 or 0.68 < L/D < 
1.25) the Cp on both upstream and downstream walls approaches -0.1, except 
near the freestream edge of the rear wall. The Cp on the rear wall in the region 


from 0.1 D to the freestream surface increases to a positive value equal to a few 
tenths. This small region accounts for the entire pressure-drag force. For cavities 
in the wake mode, shear-layer impingement occurs further below the edge of the 
rear wall. As a result the positive pressure coefficient extends much further into 
the cavity. For large L/D, impingement even occurs on the cavity floor leading to 
positive Cp values. On the upstream wall, the negative Cp is nearly twice that 


for self-sustaining oscillation modes. 


Examples of selected Cy values from Gharib and Roshko [Ref. 49] are 


presented in Figure 9(b). The plot shows the variation of maximum pressure 
coefficient with increasing cavity length. The Cp on the upstream corner of the 


cavity is nearly constant ( ~ -0.07) in the self-sustained oscillation mode II and III 


regimes and are more negative ( ~ -0.10 to -0.15) in the wake regime. In the self- 
sustaining oscillation regimes, the maximum Cp occurs at the downstream corner 


as shown by the overlap of the plus and solid-circle symbols. As cavity length 
increases in the wake-mode flow regime, the peak Cp On the downstream moves 


into the cavity from the downstream corner. This is consistent with shear-layer 
impingement inside the cavity, as discussed earlier. 
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Gharib and Roshko [Ref. 49] used a laser velocimeter to measure both the 
turbulent and mean shear layer velocity components. These proved to be difficult 
to obtain especially in the wake mode. The integral of UV across the cavity 


opening represents a relatively small difference between two large contributions. 
It was found that the variation of Cp yy, while showing a lot of scatter, was 


generally consistent with Cp p. In conclusion, it was found that the oscillating 
shear layer in the cavity is fundamentally different from a free shear layer. This 
difference is attributable to coupling between the upstream and downstream edges 


of the cavity. 
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(b) Variation of maximum pressure coefficient with cavity length L/0> [Ref. 49]. 


Figure 9. Concluded. 


In their study of coaxial flow over a combination of a disk and cylinder 
separated by a gap, Koenig and Roshko [Ref. 60] defined two flow regimes based 
on the magnitude of the drag of the combined body. These regimes were 
identified as low and medium drag. The medium drag is significantly lower than 
the drag that exists when the bodies are well separated and the downstream body 
has little influence on the upstream body, but in the low-drag regime drag is about 
an order of magnitude lower. Their study suggests that the low-drag flows appear 
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to be more general examples of the classical cavity flow. From the Gharib and 
Roshko [Ref. 49] work, it is now clear that the low-drag regime of Ref. 60 
- corresponds to the regime for self-sustained cavity flow oscillations while the 
medium-drag regime corresponds to the cavity wake mode. The flow 
visualizations and flowfield measurements of Ref. 49 show that in the low-drag 
regime the cavity shear layer always stagnates at the downstream corner. Only in 
the wake mode, the increased drag case, does it stagnate inside the downstream 
corner. The stagnation location moves further into the cavity as the cavity length 
increases. - | 


In summary, the mechanism for a stable equilibrium of the shear layer may 
be explained. An inward displacement of the shear layer.at the rear corner lowers 
the feedback ‘signal and the resultant Reynolds stress. This decreases the 
entrainment in the initial portion of the shear layer and reduces the positive 
outflow. Correspondingly, the inflow at the downstream end must be reduced and 
the shear layer returns to its initial position. In the self-sustained flow oscillation 
modes the flow is very stable and resistant to moderate external disturbances until 
the gap becomes so large that the wake-mode instability takes over. These 
phenomena have been shown to determine that a cavity in the self-sustained flow 
oscillation modes has low drag and that a cavity in the wake mode has high drag. 


D. NUMERICAL PREDICTIONS 


The earliest separated-eddy and cavity-flow computations were reviewed in 
- detail by Burggraf (Ref. 61]. These computations ranged from incompressible, 
inviscid such as Rossow [Ref. 62] to solutions of full Navier-Stokes equations at 
low Reynolds numbers. The Burggraf paper included a numerical solution to the 
full Navier-Stokes equations using a stream function and vorticity formulation, for 
the case of a square cavity in the Reynolds number (UL/v) range from 0 to 400. 
From all of the studies reviewed, a fairly clear description was obtained for the 
overall flow characteristics. However, the detailed structure of secondary vortices 
was poorly represented, even for these low Reynolds number flows. Pan and 
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Acrivos [Ref. 63] calculated steady flow in rectangular cavities (L/D from 0.25 to 
5) where the flow was driven by a uniform translation of the top wall of the cavity. 
These creeping flow solutions were obtained using a relaxation technique in a 
formulation similar to that of Burggraf. Computed streamline contours of the 
secondary corner vortices were an improvement upon earlier solutions. However, 
a minimum grid size of 0.01L was too coarse to represent detailed streamline 
patterns inside the corner eddies. It should be noted that these corner vortices 
occupied only about 0.5 percent of the total area of the cavity cross-sectional area. 
Unfortunately, the numerical method had serious instabilities for Reynolds 
numbers greater than 400, which made those results invalid. 


The flow in a two-dimensional channel with a rectangular cavity was 
studied numerically by Metha and Lavan [Ref. 64]. For simplification, the length 
of the channel was taken to be infinite and the upper wall of the channel was 
moved with a constant velocity. This kept the flow approaching the cavity 
identical in all cases. Steady laminar, incompressible flow in two-dimensional 
channels with a rectangular cutout cavity were computed using an explicit 
numerical solution of the Navier-Stokes equations in a stream function and 
vorticity formulation. Solutions were obtained for cavity L/D of 0.5, 1.0, and 2.0 
and for Reynolds numbers of 1, 10, 100, and 500. These cavity L/D ratios were 
chosen to give reattachment of the flow over the cavity and not inside. One vortex 
was observed for cavities with L/D of 1 and 2. For a L/D of 0.5 there were two 
vortices present, one on top of the other. The streamline dividing the external flow 
and the cavity flow was concave for very low Reynolds numbers and convex for 
the higher Reynolds numbers. As the Reynolds was increased, the vortex center 
moved downstream and upward, creating a thin shear layer. The shear layer on 
top of the cavity and along the cavity wall was not very thin, even at the larger 
Reynolds number. This result suggested that a relatively coarse grid may be 
adequate to define the flow features. | 


Borland [Ref. 65] obtained numerical solutions for the oscillating flowfield 
in an open cavity exposed to a high subsonic freestream for comparison with both 
wind tunnel and flight test data. Two-dimensional Euler equations for time- 
dependent inviscid compressible flow were formulated and solved using two 
different numerical algorithms. A first-order, fluid-in-cell method (FLIC) which 
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used an entirely Eulerian scheme provided the fundamental-frequency pressure- 
oscillation mode results. The grid cells were control volumes where the equations 
for conservation of mass, momentum, and energy were solved subject to 
appropriate boundary conditions. A second-order, fluid-in-cell MacCormack 
predictor-corrector method was also used in an attempt to predict some of the 
higher-frequency oscillation modes. These finite-difference equations were used 
for cell boundaries which were within the computational region. Upstream or 
inflow boundary conditions (b.c.'s) used the freestream variables. Downstream or 
outflow b.c.'s assumed that the gradients of the basic flow variables vanish. This 
is a more complex b.c.'s and sometimes produces less reliable results. The solid 
wall b.c.'s used cell boundaries which were coincident with the physical 
boundaries. For inviscid flow, reflective conditions were used where the normal 
velocities at the solid walls are cancelled by assuming a fictitious grid cell with an 
opposing velocity just beyond the physical boundary. All of the other primitive 
flow variables are also duplicated by the fictitious grid cell. The first-order FLIC 
results gave a good correlation with the fundamental oscillation frequency and the 
magnitude of the pressure oscillations observed experimentally. The second-order 
‘ predictor-corrector method showed a capability to predict higher-order harmonics 
as well as the fundamental. It was found that this computation was not stable and 
the oscillatory flow field damped out with increasing numbers of iterations and 
was not observed late in the computation. The outflow b.c. was suspected to be 
the cause of this computational instability. 


Hankey and Shang [Ref. 66] analyzed pressure oscillations in an open 
cavity using time-dependent numerical computations of supersonic flow over an 
open cavity with a L/D equal to 2.25. The unsteady Navier-Stokes equations were 
solved by the MacCormack finite-difference explicit method for a freestream 
Mach number of 1.50, and a Reynolds number of 2.6 x 10/7. The most obvious 
feature was that the flowfield within the cavity is subsonic, except in the vicinity 
of the cavity opening. The orderly development of the shear layer above the 
cavity was shown. An attached turbulent boundary layer upstream of the cavity 
separates at the front wall to form the shear layer over the cavity. The shear layer 
then reattaches downstream of the cavity. The pressure oscillation was also 
predicted, and both the fundamental frequency and magnitude were in agreement 
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with experimental data [Ref. 51]. However, the numerical solution would be 
needed for a longer period to completely verify the experimental spectral analysis. 
This paper provided the first complete viscous solution of the pressure oscillating 
cavity. 


Ghia, Ghia, and Shin [Ref. 67] used a multigrid method to solve the Navier-Stokes 
equations for incompressible flow in a cavity. The laminar incompressible flow 
was computed in a square cavity whose top wall moved with a uniform velocity in 
its own plane for moderately high values of Reynolds number (100 to 10,000). 
The objective of the paper was to achieve the solutions in a greatly reduced 
computational time. The two-dimensional cavity flow was represented 
mathematically by a stream function and vorticity form of the Navier-Stokes 
equations. The non-slip boundary conditions were applied at the cavity walls by 
requiring zero normal velocities at all of the boundaries. The two-dimensional 
cavity space was discretized by a uniform mesh. Second-order accurate central 
finite-difference approximations were employed for all of the second-order 
derivatives. The convective terms were represented by a first-order upwind 
difference scheme. In the multigrid method, the mesh density is cycled from low 
to high, iteratively, to smooth the results by damping the high frequency errors. 
The solutions used the coupled, strongly-implicit procedure of Rubin and Khosla 
[Ref. 68]. 


An example of the computational results from Reference 67 is presented in 
Figure 33 which shows the streamline contours in Figure 33 for a Reynolds 
number of 10,000 using a uniform grid (257 x 257). The center of the primary 
vortex was offset towards the top right corner at a Reynolds number of 100 (x = 
0.6172 and y = 0.7344). It moves toward the geometric center of the cavity with 
increase of the Reynolds number. Above a Reynolds number of 5000 the 
movement ceases and the primary vortex was located at x = 0.5117 and y = 
0.5233. As the Reynolds number increased from 100 to 10,000 the vorticity 
strength at the vortex center decreased from 3.17 to 1.88. 


This solution shows the presence of additional counter-rotating vortices in 


or near the cavity corners. The notation shown in Figure 10 uses the letters T, B, 
L, and R to denote top, bottom, left, and right, respectively. For example, BR 
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_ refers to the second in the sequence of secondary vortices that occur in the bottom 
right corner of the cavity. At a Reynolds number of 100 there were only two 
secondary vortices; as the Reynolds number was incresed to 10,000 the number of 
secondary vortices increased to the six shown in Figure 10. Initially, all of the 
secondary vortices are located very near the cavity corners. With increasing 
Reynolds number the secondary vortices move, very slowly, toward the cavity © 
center. Figure 33 includes magnified views of the contours of the secondary 

vortices. The values of the stream function for each contour on the figure are 


listed below: 


-0.1175 
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Figure 10. Streamline pattern for primary, secondary, and additional corner 
vortices for a Reynolds number of 10,000 using an uniform grid (257 x 257). 


In the last ten years there have been many additional computational 
solutions of the cavity flow problem [Ref. 69 through 78]. These efforts covered a 
wide range of flow conditions. In spite of this large amount of computational 
activity there are still some major deficiencies. The extensive regime of unsteady 
and oscillatory flow over and inside cavities still contains many unresolved 
problems. Prediction of the possible modes, amplitudes, and frequencies have 
been very difficult to obtain accurately. Many of these prediction efforts have not 
adequately modeled the effect of the initial boundary layer thickness on the shear 
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layer interaction with the cavity flow oscillations. Attempts to get around the 
details of the shear layer have usually succeeded only in situations where the shear 
layer is not the controlling factor, such as for deep cavities and large enclosures 
with small openings. Improvements upon the these deficiencies will require a 
significant development effort. 


E.. MULTIPLE CAVITIES | 


A two-dimensional inviscid, incompressible computational investigation of 
the aerodynamics of an airfoil with a vortex trapped by two spanwise fences was 
conducted by Rossow [Ref. 62]. To evaluate the concept, computed results were 
obtained for a Clark-Y airfoil using a variety of fence geometries at several angles- 
of-attack to determine the aerodynamic characteristics of these trapped-vortex, 
high-lift devices. These computational results suggested that two spanwise fences 
Should be used to enclose a trapped vortex. The two fence heights were adjusted 
so that an equilibrium condition could more easily be achieved than using only a 
single fence. Little or no mass removal from the core region of the vortex was 
required to stabilize the trapped vortex. The lift increase was found to be inversely — 
proportional to the chordwise spacing of the fences. It appeared that for two 
fences there would be no profile drag penalty, that the vortex would be easier to 
form and would be more stable than for a single fence. The results also showed 
that the vortex bubble could be moved fore. and aft on the airfoil to control the 
pitching moment. It was suggested that an extension of the single trapped-vortex 
geometry to two or more trapped vortices on the upper surface of an airfoil may 
provide the same lift with less cumbersome. equipment. 


The progression in complexity from the two-dimensional configurations 
studied above to the full three-dimensional flowfield of a wing requires 
examination of several additional factors. Such devices were recently studied in 
an unpublished experimental investigation by Rossow and Ross [Ref. 79] using a 
low-aspect-ratio, highly-swept-back (A = 60°) wing (Figure 11) with a semispan 
of 24.76 inches and a chord of 18.98 inches. There were differences between the 
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actual wing area and the resultant aspect ratio and the values used in the 
aerodynamic coefficients as the reference geometry as shown below: 


a ee 


This test was conducted in the NASA Ames Research Center 7- by 10- Foot 
Tunnel using a semispan model. The model was mounted on a metric balance 
frame. Force data were obtained from an external floating ame system of beam 
balances located below the test-section floor. 







V 


———— 


Figure 11. Semispan wing (NACA 0012 airfoil normal to the leading edge) with 


two part-span spoilers tested by Rossow and Ross [Ref. 79]. 


A variety of spoiler configurations were tested. As a sample of results from 
this investigation, the following data were obtained using a pair of spoilers (0.08c 
high at 0.15c and 0.10c high at 0.45c) which were used to form a single cavity 
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swept along the local chord line. The spoilers extended from the root chord to 
approximately 1/3 of the wing span. The intent was to take advantage of the 
inherent spanwise flow to form a standing vortex and to avoid active mass removal 
from the core region of the vortex. If the spoilers are effective, then the flow over 
the upper surface would be displaced passively to increase the wing camber. This 
experimental investigation found that below 8° angle-of-attack there was little or 
no lift increase. However, there was a lift benefit as high as 30% at an angle-of- 
attack greater than 10° (Figure 12(a)). 





a, deg. 


(a) Lift curve. 
Figure 12. Rossow-Ross experimental data [Ref. 79]. 


However, there was an unexpected and undesired drag increase over a wider 
angle-of-attack range as shown by the polar plot in Figure 12(b). For reference, 
the ideal polar (1/m AR e) is shown by a solid line, where the span efficiency . 
factor e is 1, which corresponds to an elliptic spanload distribution. An efficient 
transport wing will typically have a span efficiency factor greater than 0.9. 
Lower-aspect-ratio fighter aircraft have a span efficiency which drops into the 0.7 
or 0.8 range. Even lower span-efficiency factors indicate either an inefficient 
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spanload distribution or the presence of flow separation.The polar for the wing 
without spoilers shows a drag increase above the ideal drag polar for Cy, > 0.3 (or 


o> 5°). The drag for the spoiler-deployed configuration is increased even further. 
In an attempt to quantify the drag increase the experimentally determined ven 
efficiency factor, e, is tabulated in the following table: 


ig 


seen fl ; cua uae a wollen deployed 
3 : none - ideal polar 
| — O-— none : 
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These adverse drag results suggest the presence of flow separation and indicate the 
need for a better understanding of the flow phenomena involved. 


In the early 1960's, Migay [Ref. 80] conducted a series of investigations 
which evaluated the effect of multiple transverse cavities in a diffuser (Figure 
13(a)). The diffuser surface had a constant radius of 11 centimeters. The first 
cavity was located at the entrance to the diffuser ( @ = 0°). Up to 12 transverse 
cavities distributed along the flow direction through the diffuser (at « > 0°) were 
tested. Pressure taps were distributed along the diffuser in the direction of the 
flow. The pressures were non-dimensionalized using a pressure coefficient 
referenced to atmospheric pressure : 


e= Pi Patm (12) 
ps] 2 
5 (Ule-0) 


where U|,_, is the velocity in the test section upstream of the diffuser. The 


pressure data were presented as the ratio of the pressure coefficient measured at a 
particular tap location to the pressure coefficient measured at & = 0. The measured 


pressure distribution is presented in Figure 13 (b) and it is seen that the largest 


curved diffuser with 

15 cavities located 3.5 mm apart 
each cavity is - 

7 mm deep & 1.5 mm wide 


70 mm 


110 mm radius 


(a) Sketch of curved diffuser. 
Figure 13. Results of Migay investigation [Ref. 80]. 
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(b) Effect of diffuser angular position. 
Figure 13. Continued. 
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(c) Effect of the total number of cavities 
Figure 13. Concluded. 
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pressure reduction occured near a = 10°. It appeared that the surface flow 
separated downstream of the & = 10° location. It was found that the cavities were 
effective in preventing flow separation as long as the first cavity was located at an 
a between 0° and 10°. When the first cavity was located further downstream, the 
cavities had no beneficial effect on the location of flow separation. The effect of 
the number of cavities is presented in Figure 13 (c) where the ratio of pressure 
coefficients is shown as a function of the number of cavities. No effect was found 
for up to three cavities. Then there was an increasing pressure coefficient 
reduction from four to eight cavities with no additional change for more than eight 
cavities. These results indicate that the cavity benefit can be achieved with a 
limited number of cavities if they are properly located in the diffuser. 


An investigation by Lin et al [Ref. 81-86] evaluated in some detail the 
potential of multiple cavities as a base-drag reduction device. An example of these 
results is presented in Figure 14 for an axisymmetric body with a modified bluff 
base which had up to 13 circumferential grooves. The transition from the 
axisymmetric body to the tapered base used a radius on the shoulder at the 
intersection of the body and the base. For the results in Figure 14, the shoulder 
radius was one-half of the body radius. The downstream edge of the first groove 
was located at the start of the shoulder radius. The grooves were numbered from 1 
to 13 as their location varied from the shoulder to more downstream positions. 
The bluff-body base had a 30° angle; a streamline body with a base angle of 10° 
was also tested. Both bodies had a diameter of 2.39 inches; the bluff body was 3 
diameters long while the streamline body was nearly 5 diameters long. 


The data in Figure 14 present the variation of drag coefficient with 
Reynolds number for axisymmetric-body configurations with and without 
circumferential grooves. They show that the streamline body had the lowest drag 
and that the ungrooved bluff body had one of the highest drag coefficients. The 
effect of the grooves was determined by opening the grooves one at a time. 
Groove number 1 was opened at the shoulder location. Then groove number 2 
was also opened at the next downstream location. This process was continued 
_ until all 13 grooves were open. When one groove was opened the drag increased 
above that for the ungrooved body. When a second groove was opened there was 
a drag reduction. However, for 1, 2, or 3 grooves open, the drag was greater than 
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for the ungrooved body. As the number of grooves was increased from 4 to 10 the 
drag was increasingly reduced below that for the ungrooved body. 


The largest drag reduction was achieved with either 10 or 13 grooves open. 
For these two configurations the lowest drag occurred in the Reynolds number 
range from about 80,000 to 100,000 and then drag tended to increase slowly with 
further Reynolds number increase. It was found that the boundary layer 
transitioned on the base from laminar to turbulent in this Reynolds number range. 
As a result, these circumferential grooves on the body with a shoulder radius 
appear to only be effective for a laminar boundary layer approaching the shoulder. 
The reason for this Reynolds-number sensitivity was not identified. This 
investigation also evaluated many other devices and found several which were as 
good or better than multiple cavities without being limited to a particular 
Reynolds-number range. This study also indicated a need to improve the 
understanding of the flow changes due to the cavities. 
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Figure 14. Effect of circumferencial grooves on the afterbody on the 
axisymmetric bluff-body drag; shoulder radius is 0.2 body radius [Ref. 86]. 
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il. EXPERIMENTAL PROCEDURE AND RESULTS 


The present investigation included a wind tunnel test where the effect of 
single and multiple cavities on the boundary layer development and the pressure 
distributions were measured in both a zero and two adverse pressure gradients. 
The boundary layer profiles were measured at four or five longitudinal stations in 
the vicinity of the cavity(s). One station was located immediately upstream of the 
first cavity and the rest were either between the multiple cavities and/or 
downstream of the cavities. 


A. WIND TUNNEL AND APPARATUS 


The experimental investigation was conducted in the San Jose State 
University (SJSU) Low Speed Wind Tunnel. The test apparatus included a three- 
dimensional traverse from the Naval Postgraduate School as well as a computer 
and data aquisition system from the NASA Ames Research Center. The test 
hardware for the cavity configurations and the test section diffusers were 

specifically fabricated for this investigation in the SJSU tunnel. | 


1. Wind Tunnel 


A sketch of the SJSU tunnel is presented in Figure 15. This is a continuous 
flow tunnel with an Eiffel-type return. The flow enters the 30" square entrance 
section, shown on the left in Figure 16, and then goes through a honeycomb and 
two screens into the 36" long contraction. The contraction ratio is 6.25. The test 
section is 12" high, 12" wide, and 24" long. The contraction, test section, and 
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diffuser are mounted on a fixed frame. The drive motor-fan assembly is located 
on a movable frame located on the right side of the sketch. There is also a 
honeycomb upstream of the motor. Between the downstream end of the diffuser 
and the honeycomb there is a variable gap which is used to control the tunnel 
speed. As the gap was varied from 6" to 0" the test section speed Uso increased 
from 40 feet per second (fps) (q = 1.9 psf) to the maximum tunnel speed, 142 fps 
_(q = 24 psf). The present test was conducted using a speed of 130 feet per second. 
A 96" long test section extension was constructed for this investigation to increase 
the boundary layer thickness for a portion of the tests. Longitudinal position in the 
test section was measured from the test section entrance: station 0" through 24" for 
the 24" test section; and stations 0" through 120" for the 120" test section. In the 
longer test section the original 24" long test section was located between the 96" 
and 120" stations. Photographs of the SJSU tunnel with each test section length 
are presented in Figure 17. 


To achieve adverse pressure gradients in the 24" long test section, the first 
6" was reduced in height. Then the increased contraction ratio was faired using an 
insert on the tunnel bottom in the 16" immediately upstream of the test section 
(x<O0"). A test section diffuser was installed between 6" and 24" stations on the 
floor of the test section. One diffuser had a constant 7° angle and the second test 
section diffuser had a constant 14° angle. The test section extension height was 
variable to accomodate the reduced height upstream of the same test section 
diffusers used in the 24" test section. The contraction ratios for the test section 
with the 7° and 14° test section diffusers were 9.38 and 10.13, respectively. 
Coordinates are presented in Figure 36 for the original test section (solid lines) and 
for the two test section diffusers (dashed lines). The x and y scales exagerate the y 
dimensions relative to the x dimensions. Tests were conducted in each of the three 
configurations shown; (1) original test section, y = 6" (top) and -6"( bottom); (2) | 
7° test section diffuser, y = 6" (top) and varying from -4" to -6" (bottom); and 14° 
test section diffuser, y = 6" (top) and varying from -1.4" to -6" (bottom). The 
same test section top was used for each of the basic tunnel configurations. For the 
120" test section the same test section diffusers were used. As mentioned in the 
previous paragraph their x locations were increased by 96", the length of the test — 
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Figure 15. Perspective sketch of theSan Jose State University Low Speed 
Tunnel. : 
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(a) Coordinates of the 24 " test section and for the 7° and 14° test section diffusers. 
Figure 18. Test section diffusers used to generate an adverse dp/dx. 


section extension which was added to the front of the original test section and 
downstream of the modified contraction ratio (x < 0"). 


Photographs of the SJSU tunnel are also presented in Figure 18 which 
shows the 7° test section diffuser with four one inch cavities (Figure 18(b)) and the 
14° test section diffuser with one two inch cavity (Figure 18(c)). The test section 
diffusers were fabricated using plywood frames mounted in the streamwise 
direction and covered by 0.25" thick clear acrylic sheets. The photograph for the 
7° test section diffuser also shows the four cavity configuration mounted on the 
top of the test section. The photograph for the 14° test section diffuser also shows 
the single two inch cavity configuration mounted on the top of the test section. 
Some of the clutter in the background of the photographs is the tygon tubing used 
for the pressure distribution measurements on the tunnel top and bottom. The 
vertical portion of the pitot-static probe used for tunnel reference dynamic pressure 
measurement is visible above the circular access window. The boundary layer 
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total pressure probe is near the downstream cavity in Figure 18(b) or near the two 
inch cavity in Figure 18(c) adjacent to the tunnel ceiling in both photographs. 


2. Test Hardware 


The test hardware built assembled specifically for the present investigation 
consisted of the test section diffuser, the cavity configurations (1, 2, and 4 
cavities), the three-dimensional traverse mechanism, and the computer with its 
data acquisition program. | 


a. Single and Multiple Cavity Hardware 


There were four cavity configurations built for the investigation of 
the effect of the number of cavities on the boundary layer profiles. As shown by 
the photographs in Figure 18 the cavities were mounted in the test section ceiling 
in the adverse pressure gradient portion of the test section diffuser. For 
nomenclature purposes, the edge of the cavity adjacent to the freestream is called 
the top and the cavity horizontal surface parallel to the freestream is called. the 
floor. For each cavity the upstream vertical surface is called the front wall and the 
downstream vertical surface is called the rear wall. Each cavity had a nominal 
length to depth ratio L/D of 1. Each of these cavity configurations was equipped 
with a distribution of surface pressure taps both in the cavities and on the ceiling 
surface both upstream and downstream of the cavities. The coordinates of all of 
these pressure taps are presented in Appendix B, Pressure Tap Locations. The 
cavity test hardware which was fabricated for this investigation is summarized in 
the following table: | 
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number of | cavity number cavity length, | longitudinal 


cavity depth, 
location of 





cavities (from 


upstream) 





cavity front 
wall, in. 





11.95 


b. Traverse 


The three-dimensional traverse assembly from the Naval 
Postgraduate School is shown in Figure 19: It was used to support and traverse a 
pitot probe used to measure the boundary layer profile. A 0.375 inch circular hole 
in the ceiling provided clearance for the 0.25 inch diameter probe support. The 
probe was used to survey the boundary layer profile adjacent to the test section 
ceiling. The traverse was mounted on top of a frame which wrapped around the 
tunnel test section and did not touch the wind tunnel. Three electric motors 
powered the gear-driven screws which moved the traverse independently in each 
Cartesian direction. The following traverse ranges were available: streamwise, 
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Figure 19. Photograph of the traverse mechanism mounted at the SJSU tunnel. 


0"; laterally, 22.5"; and vertically, 12". These movements were executed by 
manually entering the desired movement direction (X, Y, or Z) and the movement 
distance. The repeatibility of probe positioning was determined by moving the 
probe ten inches and then moving back to the starting point. There was less than 
0.002 inch error in returning to the original position. Significantly no backlash 
due to reversing the movement direction was observed. To assure the best 
positioning accuracy, all of the boundary layer traverses were obtained by moving 
in one direction from the ceiling into the freestream to avoid any possibility of 
backlash. For the next run, the probe was then returned to a position hard against 





the test section ceiling. The direction was then reversed to a position where the 
probe just touched the ceiling prior to the next traverse. 


c. Boundary Layer Probe 


The boundary layer total pressure was measured using a total 
pressure probe which was made by flattening a circular stainless steel tube. Using 
a 0.060 inch diameter tube, a probe tip height of 0.022 inch was used for runs 101 
through 254. This probe was damaged and replaced by another probe with a tip 
height of 0.034 inch for runs 255 through 380. The probe location for plotting 
boundary layer profiles was taken to be one-half of the probe tip height when it 
touched the wall. For the tests, this half height was identified as zero on the probe 
traverse position data display. 


The pressure was measured using a 0.5 psid transducer. The NASA 
Ames calibration laboratory found that this particular transducer had a standard 
deviation equal to 0.141% of full scale. This translated to a 0.5% (0.102 psf) error 
for differential pressures measured at a freestream dynamic pressure of 20 psf. 
After the probe was installed, the calibration was checked using an electronic 
pressure calibration device. Regression analyzes of repeated calibrations gave a 
correlation coefficient r2 of between 0.9986 and 0.9997. This represented a 
maximum repeatibility error of 0.07% of full scale or 0.05 psf. A third error 
source was the electronic zero shift between the beginning and end of arun. In 
154 runs this shift averaged 0.03 psf. If these errors were summed the maximum 
error becomes 0.18 psf or 0.9% of the freestream dynamic pressure. 


A final error source was the effect of the wall proximity or boundary 
layer displacement effect. The work of Coles [Ref. 9] was applied to the present 
data. It was found that this error was negligible when used to correct the 
displacement thickness 6] and the momentum thickness 62. Application of this 
correction to the measured velocity adjacent to the wall indicated a velocity error 
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_AV/V of 0.0145. At a freestream dynamic pressure of 20 psf, this was an error of 
0.004 psf or 0.02%. This was considered to be negligible. 


d. Computer and Data Acquisition System 


The digital data acquisition system was controlled by a desktop 
computer through a graphical user interface (GUI) computer program and two 
dedicated data acquisition circuit boards. Each data channel had its own amplifier 

- which conditioned the data signal and amplified its voltage so that it could be 
sampled with the greatest possible accuracy and then multiplexed to the computer. 
The computer had 8 megabytes (Mb) of random access memory (RAM), an 40 Mb © 
internal hard disk, and a removable 44 Mb hard disk which was used to store the 
experimental data. 


The application program for the GUI was a set of virtual instruments 
(VI) which each consistd of three parts: (1) the front panel; (2) the block diagram; 
and (3) the icon/connector. The first two parts were visible on the computer 
screen as windows. The front panel was the user interface to the VI and was used 
during data acquisition by the program operator to interact with it. It included 
input objects called controls and output objects called indicators. The objects were 
portrayed on the computer screen as command buttons and status indicators. Data 
was also be displayed in either tabular and/or grapical forms. The block diagram 
was a VI source code which was created using a graphical programming language. 
This was a virtual wiring diagram for connecting the program modules which were 
called nodes to the terminals for the controls and indicators. The icon was a 
representation of the VI. The connector actd as the port through which the data 
passed. The icon/connector together represented the VI in a manner analogous to 
a subroutine call statement when the VI was used as a subVI in another VI's block 
diagram. Many layers of subVIs could be used in higher order VIs. 


The present investigation GUI used a main panel (Figure 20(a)), a 
boundary layer measurements panel (Figure 20(b)), and a pressure distribution 


56 














_ measurements panel (Figure 20(c)). The data acquisition/reduction program was 
executed by first loading the GUI program application code and then launching the 
cavity flow experiment main panel VI and all of the related subVIs. At the 
completion of these steps the computer screen displayed the main panel shown in 
Figure 20(a). There were nine command buttons shown on the top portion of this 
panel and status indicators for each data channel. The data acquisition was 
initiated by sequentially using the computer's mouse to select the clear acquisition 
button, then select the initialize acquisiton buttom, and finally select the take zeros 
button. Then there were five choices available: (1) data acquisition setup button 
which was used to add or remove data channels; (2) on-line monitoring button to 
view data output for each channel; (3) boundary layer measurements button; (4) 
pressure distribution measurements button; and (5) quit button to terminate and 
exit the program. | 


In this investigation four groups of data measurements were taken 
and recorded in data files: (a) tunnel conditions; (b) boundary layer profiles; and 
(c) static pressure distributions. An example of the data recorded for a boundary 
layer profile run is presented in Figure 21. This particular run was identified as 
RO1OSBL which stands for run R, number 0105 and boundary layer profile BL, 
respectively. Subsequent boundary layer profile runs have different run numbers. 
As shown in Figure 21, the data file generated by the GUI program started with the 
run number, the date and time of the run, and the tunnel code (see Appendix A for 
a detailed description). Then there was an additional line of group 1 data which 
consisted of the freestream conditions measured by a pitot probe located at x = 1" 
and y = 2.5" and then averaged over the entire run. These data include the tunnel 
reference values for pitot probe static Ps,tun and total pressures Pt,tun , as well 
as their difference, the calculated dynamic pressure gc. For an independent 
measurement, the dynamic pressure was also measured by a separate pressure 
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transducer and displayed as q. For pressure coefficients this value of dynamic 
pressure was used to non-dimensionalize the static pressures. Finally the tunnel 
total temperature Tt,tun and atmospheric pressure P atm were listed. The tunnel 
conditions were always measured along with one of the other three groups of data 
described above. 


Shown in Figure 21 is a group of boundary layer probe data, 
identified as group 2 data. The number of data points varied from run to run. In 
this example from run 105, 15 points were recorded. Then seven columns of data 
measurements were displayed for each boundary layer height: (1) the height from 
the bottom of the boundary layer total pressure probe to the tunnel ceiling (y = 
6.00") Y bl; (2) the boundary layer total pressure Pt,p; (3) the ratio of local u 
velocity to estimated local freestream velocity u/U; (4) pitot probe dynamic 
pressure q,tun; (5) tunnel total temperature Tt,tun; (6) pitot probe total pressure 
Pt,tun; and (7) pitot probe static pressure Ps,p. The height of the boundary layer 
measurement was the sum of the probe height Y bl plus one-half of the probe 
vertical thickness Ay)/2 which was 0.011 inch for runs 100-257 and 0.017 inch 
for runs 258-380. The four measured pressures and the measured temperature 
were tabulated in the engineering units shown in Figure 21. The boundary layer 
profile in run 105 was taken at x = 9", which is downstream of the location of the 
single cavity configurations. The boundary layer profile is expressed as the ratio 
of boundary layer velocity computed from the total pressure measured by the 
boundary layer probe Ptp at the height Y bl + AYp/2 divided by the freestream 
velocity at the x station where the profile was measured. The velocity within the 
boundary layer at a given height was computed from the Bernoulli equation as 


2 33 
u=. 5 (Pep ~ Ps,ceiling J (13) 


where the static pressure measured at a static pressure orifice located on the ceiling 
at the x station of the boundary layer profile was poy,¢. The quantity Pt.p-Psurf 


was measured as the pressure difference across a pressure transducer. 


During the on-line measurements the local freestream velocity was 
estimated by the GUI program. For tunnel configurations without the diffuser 
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inserts, the local freestream velocity was assumed to be constant along the length 
of the test section. For tunnel configurations with a diffuser the local velocity U 
or Ux was computed using the incompressible continuity equation 
U, =U, Aco (14) 
A, 
where either; (1) A, =12*y@,-5 for A,, or for Ay where x < 6 or; (2) 
A, = 12(y @xa5 t(X— 6) sin(o.) for A, when x > 6. The test section diffuser had 


a constant height from x = 0" to x = 6". From x = 6" to x = 24" the diffuser had a 
constant slope of either 7° or 14°. For the final processed data, Ux was replaced 
by the velocity at the edge of the boundary layer Up. which was computed from 
the boundary layer probe measurements at y 2 d where the maximum value of Pt.p 
- Poyrf Was measured. 


An example of group 3 data from run 173 is presented in Figure 22. 
This particular run was identified as RO173CP which stands for run R, run 
number 0173 and pressure coefficient data CP, respectively. The run heading, 
run-date-time-tunnel code, and group 1 data are displayed the same way they were 
shown in Figure 21. The group 3 data, pressure distributions, are grouped by the 
scanivalve module that measured a group of pressure taps. Each section of these 
data is headed by the scanivalve number (1, 2, or 3) and the number of pressures 
recorded. Notice that the number of pressures is one greater than the number of 
taps for each scanivalve. The extra pressure is a reference pressure which 
represents zero return after the group of taps were sampled and indicates whether 
there is any electrical drift in the pressure data during the data sample period. 


Scanivalve | (Figure 22(a)) measured the 16 taps located on the test 
section diffuser between x =~ 2" and x = 21". When there was no test section 
diffuser installed, scanivalve 1 was not sampled; only the test section ceiling 
pressures were measured. For run 173 scanivalve 2 (Figure 22(b)) measured the 
17 taps located on the test section ceiling upstream (from x = 2") and downstream 
of the two inch cavity (to x ~ 22"). Notice that pressure tap 13 measured a 
pressure nearly equal to zero for this run. Since this was not seen in measurements 
of this tap in other runs, this measurement was interpreted as faulty data for this 
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run and omitted from the final data plots. Scanivalve 3 (Figure 22(c)) measured 
the 21 taps located in and near the two inch cavity. As seen in the data listing, the 
upstream cavity wall was located at x = 7" and the downstream wall was located at 
x = 9". The floor of the cavity was located at z = 8" (for the 2" cavity). The data 
listing for each scanivalve includes the following information: (1) tap number; (2) 
tap x location; (3) tap y location; (4) the surface static pressure Psurf; (5) the 
pressure coefficient Cp which is defined as 


C, _ Psurf ~ Ps,tun | (15) 
q 


The complete list of runs completed in this investigation is presented in Appendix 
A. The surface tap locations for all of the test configurations are listed in 
Appendix B. The cavity pressure data are presented in Appendix C for all of the 
measured data. The boundary layer profile data are presented in Appendix D for 
all of the measured data. | 
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B. WIND TUNNEL RESULTS 


_ The experimental wind tunnel data described in the previous section and 
presented in the appendices are analyzed in this section in three parts: (1) 
boundary layer profile analyses; (2) static surface and cavity pressure data 
analyses; and (3) dynamic flow due to cavity opening. The entrance velocity 
profile was surveyed in the 24 inch test section at the x = 4" station using the 
boundary layer total pressure probe. It was found that the flow dynamic pressure 
was uniform (less than 0.2 percent variation) with an increase of about 0.8 percent 
in dynamic pressure at the edge of the ceiling and floor boundary layers. A single 
hot-wire probe was located at the tunnel centerline of the 24" test section at the x = 
4" station to measure the streamwise flow turbulence. These data found that the 
Au/u was approximately 0.003. 


1. Test Section Boundary Layer 


For the plain test section ceiling, the boundary layers were measured at four 
longitudinal locations in the original 24" test section (x = 5.5", 9.0", 12.0", and 
16") and in the extended 120" test section (x = 101.5", 105.0", 108.0", and 112.0"). 
The test-section configurations included the basic 12" by 12" cross section test 
section with essentially zero pressure gradient (dp/dx ~ 0) and the two test section 
diffuser configurations with adverse pressure gradients (dp/dx > 0). These data are 
analyzed in the following section to evaluate how representative the present basic 
tunnel data are of a two-dimensional boundary layer over a flat plate. 


a. Zero Longitudinal Pressure Gradient 


For a two-dimensional, turbulent boundary layer in a zero pressure 
gradient, the velocity profile may be estimated using the 1/7 power law 
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u/U,, =(y/ §)\" ” The 1/7 power law may then be integrated [Ref. 4] to obtain the 
above boundary layer thickness quantities as 96/x=0.37 R,/) 
6,/x = 0.37 Ro") /(1+n) where n = 7, and 8,/x=0.036 R,~!9). In the 
present investigation these relations only apply to the cases where the velocity is 


constant with increasing run length x ; that is for zero pressure gradient, the tunnel 
configurations without the 7° and 14° diffusers in the test section. | 


These three boundary layer thickness parameters (6, 61, and 69) 


were computed using both the 1/7 power law relations and the experimental data 
from the SJSU tunnel for both the basic 24" long test section (at x = 5.5", 9", 12", | 
and 16") and the 120" long test section (at x = 101.5", 105", 108", and 112"). The 
results are compared in Figure 23. In the 24" test section data all three 
experimental thicknesses were found to be slightly greater than the 1/7 power law 
thicknesses. This could be due to the fact that x was measured from the start of the 
test section. Ahead of the test section in the contraction the pressure was . 
decreasing (dp/dx < 0) so that the boundary layer thicknesses were decreasing. 
However, they almost certainly do not go to zero at the x = 0 station. If the initial 
station were moved upstream 3 or 4 inches to a virtual origin for application of the 
1/7 power law (to more adequately represent the equivalent zero thickness 
location) then the experimental data would agree with the 1/7 power law relations. 
The experimental boundary layer profiles for all four stations where the boundary 
layer was measured in the 24" test section agree with the 1/7 power law profile. 
As the Reynolds number increases the boundary layer profiles are more nearly 
one-eighth or one-ninth or a lower power profile. In the 120" test section the 
boundary layer thicknesses are well below those obtained from the 1/7 power law. 
If a one-eleventh power law is used to compute the momentum thickness 0, 
agreement with the experimental data can be achieved. For the thicker boundary 
_ layer in the 120" test section, the presence of four walls may be the reason that the 
boundary layer thickness does not grow as large as a two-dimensional boundary 
layer. While all of the reasons for the difference between the 24" and 120" test 
sections have not been identified, these effects, especially for the longer test 
section, will be evaluated numerically in a later section using the incompressible 
Navier-Stokes equations. 
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(a) Boundary-layer thickness, 6. 
Figure 23. Comparison of experimental data from present investigation using the 
24" and 120" test sections with the 1/7 power-law boundary-layer thicknesses. 
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(b) Boundary-layer-displacement thickness, 61. 
Figure 23. Continued. 
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_(c) Boundary-layer-momentum thickness, 65. 
Figure 23. Concluded. 


. Rotta's analysis [Ref. 9] introduced universal similarity parameters 
which were later demonstrated by Clauser [Ref. 8]. Clauser plotted data from 
three separate investigations as u/U versus y/d which appeared to lack similarity. 
He then replotted these data using a universal similarity-parameter plot ((u-U)/u* 
versus y/O). This change greatly reduced the data scatter and convincingly 
demonstrated similarity independent of both Reynolds number and roughness. A 
similar correlation was done using the data from the present investigation where 
the 24" test section data (Figure D1) deviated from the 1/7 power law and the 120" 
test section data (Figure D2) tended to agree with the 1/7 power law. The 
universal similarity parameter form of these data is presented in Figure 24(a). The 
data provide a good correlation for most of the boundary-layer profile with only a 
limited dispersion at heights below 0.15 for the shortest run length data (x = 5.5 
inches). This correlation provides confidence that these boundary layer profiles 
are similar and representative of those measured in a two-dimensional boundary- 
layer. 
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(a) Experimental data from present investigation. 
Figure 24. Universal velocity distribution law comparisons. 


In Schlichting [Ref. 4] the velocity defect laws of both von Karman 


u-U. My f_¥ le fh 
= Auk 1 aa i (16) 


Ueru_t ft 





and Prandtl 





compared favorably with experimental data from both smooth and rough pipes. It 
is stated that both laws were obtained for two-dimensional flow in a channel. 
They were found to agree with both two-dimensional and axisymmetric 
experimental data. Both of these equations are compared in Figure 24(b) with zero 
pressure gradient, flat plate experimental data both from the present investigation 
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(Figure 24(a)) and from the experimental data correlation done by Clauser [Ref. 
8|. The two sets of experimental data correlate. In contrast to Ref. 4, there is a 
difference between the two velocity defect laws (equations 16 and 17) and the two 
sets of experimental data in Figure 24 (b). 


: 3 Prandtl 
(U-u)/u* LY rt. 3 : ------ von Kérman 





y/d 


(b) Prandtl and von Kaérmdén equations compared with present investigation and 
with expermental data correlated by Clauser [Ref. 17]. 


Figure 24. Concluded. 


The data from the SJSU low speed tunnel is compared with log-law 
relation 


ut =<In(y*)+C (18) 


in Figure 25 (a) for the standard 24" long test section (where 6 = 0.23 at x = 5.5" 
and 6 = 0.45" at x = 16") and in Figure 25 (b) for the extended 120” long test 
section (6 ~ 1.17" at the locations where the measurements were taken). The 
comparison is typical of the expected agreement in the log-law region. There is a 
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(a) 24" test section 
Figure 25. Log-law relation using inner variables compared with experimental 
data from the San Jose State University low speed tunnel. 
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(b) 120" test section. 
Figure 25. Concluded. 
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boundary layer probe displacement effect at y+ for low heights in the inner region; 
and there is a wake effect at y+ for the greater heights. The wake is characterized 
by a region at heights approaching the edge of the boundary layer region where the 
measured velocity is greater than that given by the log law. The specific locations 
for these effects are given below: 


die oY ml 
boundary layer probe displacement effect below yt 100 | 100 
boundary layer wake effect at yt greater than 1000 






In Figure 25 (a) for the 24” test section there is a larger wake effect 
at the x = 12" location than at the other locations. It should be noted that the 
boundary layer run lengths are relatively short and so it should not be surprising 
that there are differences where the boundary layer is in its initial development. In 
contrast, as shown in Figure 25 (b) for the 120" test section, where there is a fully 
developed turbulent boundary layer, there are no significant differences among the 
wakes. As a consequence, for the y* regions both below the boundary layer 
displacement effect limit and above the wake effect limit the experimental Ut 
data is greater than the log law relation. In between these yt limits where 


equation 13 applies, the data in the log law region agrees well with the above log- 
law equation when the constants are k = 0.41 and Cy = 5.0. 


Computation of the turbulent boundary layer was the subject of the 1968 
AFSOR-IFP-Stanford Conference on Computation of Turbulent Boundary Layers 
[Ref. 9]. Volume II of reference 9 included a compilation of experimental 
boundary layer data by Coles and Hirst and a paper by Coles which was a guide to 
the experimental data compiled for use at the conference by the various predictive 
methods which were being evaluated. In this paper, Coles discussed two problems 
associated with evaluation of experimental data: (1) the boundary layer 
displacement effect where, near the wall, the data are in error because of the 
effects of probe wall interference and local strong turbulence; and (2) the wake 


P, 














effect near the edge of the boundary layer where the data are correct but the law- 
of-the-wake equations are not correct as discussed earlier in this section. To 
account for the near wall effects, Coles used a standard sublayer profile to obtain 
the following standard functions to reproduce the tabulated thicknesses: 


f. Utay* = 540.6 and j. (Ut) ay* =6546. (19) 


The experimental boundary layer displacement thickness and the momentum 
thickness have been recomputed using these relations (equations 19) and are 
shown in Figures 26(a) and 26(b), respectively. Here the 1/7 power law results are 
compared with the experimental data from the SJSU tunnel for both the integrated 
uncorrected data and these data modified by the above correction procedure. It is 
shown that the corrections made only a small change to 61 and 09. Since these 
changes are so small, for all of the rest of the boundary layer data in the present 
paper only the integrated experimental data without the Coles correction will be 
used. 


—o— 1/7 power law 
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(a) Boundary-layer displacement-thickness. 
Figure 26. Comparison of 1/7 power law boundary-layer thicknesses with 
experimental data both with and without Coles near-wall correction. 
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(b) Boundary-layer momentum-thickness. 
Figure 26. Concluded. 


Relations for turbulent wall-friction were determined using the 
experimental velocity profile data. Experimentally, the skin friction coefficient 
may be found using a variety of skin-friction strain gauge-devices, by an obstacle 
block, by a heat transfer analogy [Ref. 16], by surface flow visualization, by a 
razor blade technique, by a sublayer fence [Ref. 14 through 17], by a Preston tube 
[Ref. 18], by empirical correlation [Ref. 17 and 19] or by analysis of the boundary 
layer profile [Ref. 8, 17, and 20]. In the present report the skin friction was 
determined from the boundary layer profile using a Clauser plot [Ref. 8]. The log- 
law equation (equation 8) may be re-arranged to become an implicit equation for 
skin friction coefficient, giving 


Re +—In n==2h We lec (20) 


An example of the application of this equation to the analysis of the data from the 
present investigation is presented in Figure 27 for station 112 inches in the 120 


TT 





log (R,) 


. Figure 27. Application of the Clauser plot for the determination of skin-friction 
coefficient for data from SJSU 120" test section at station 112 inches. 


inch long test section where the skin-friction coefficient was found to be 0.0030. 
Typically, repeat runs had differences in cf no larger than 0.0002. These results 
show that while there is a small reduction in skin friction due to the cavity, the 
differences from the plain wall configuration were nearly within measurement 
accuracy. This procedure was used to obtain the skin-friction coefficient from all 
of the data described in the test run schedule (Appendix A). 


There are a number of correlations which may be used to estimate the skin- 
friction coefficient on a flat plate. The equation based on the 1/7 power law has 
already been given as cr =0.0592(R, ye . Another correlation often used is 


the empirically developed skin-friction equation by Schultz-Grunow [Ref. 19] 
ce = 0.37(logR, ) 7 : (21) 


The skin friction coefficients obtained from the Clauser plots for the conditions 
with zero pressure gradient are compared in Figure 28 (a) with the corresponding 
values obtained from the 1/7 power law and with the Schultz-Grunow equation. 
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The correlations and the experimental data agree for the 24" test section. But for 
the 120" test section both cr correlations are about 0.0005 less than the 
experimental data. Whether this difference may be due to surface roughness in the 
96" test section extension will be evaluated in the fourth subsection of the general 
background section. | 


The following correlation by Ludwieg-Tillmann [Ref. 17] was developed to 
account for the effect of pressure gradients 


c, = 0.246(Rs, i 


10)70-678H (22) 
This correlation is compared with data from both test section lengths in Figure 28 
(b). In contrast to the 1/7 power law and Schultz-Grunow correlations in Figure 
28 (a), the Ludwieg-Tillmann correlation demonstrates good agreement with the 
Clauser method results. There is scatter for the data from the 24" test section. 
There is agreement with the experimental data from the 120" test section. This 
latter comparison for x > 100" represents the major difference between the 
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(a) Comparison of Clauser analysis of present data from both the 14" and 120" 
test sections with correlations from 1/7 power law and by Schultz-Grunow. 


Figure 28. Comparisons of skin friction from several methods. 
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(b) Comparison of Clauser analysis with the Ludwieg-Tillmann correlation. 
Figure 28. Concluded. : 


Ludwieg-Tillmann correlation and both the 1/7 power law and Schultz-Grunow 
correlation. 


Initially, the pressure distributions for the 24" and 120" long test 
sections were measured. An example of the measured pressure distribution on the 


short test section ceiling is presented in Figure 29 using open circle symbols. 
These data show a reduction in Cy with increasing x distance from the test section 


dc 
entrance. The measured pressure gradient was — =—().000643. This Cp 
X 
decrease was caused by the small increase in local velocity due to the expected 
increase in displacement thickness. At an arbitrary x location the effective test 
section cross-section area for the flow may be expressed as the physical dimension 


reduced by the boundary layer displacement thickness on each surface, giving — 
A=(12.0-26,)° (23) 


and the local velocity using the incompressible continuity equation becomes 
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U, = U,A, / A» | (24) 


Finally, the computed Cp is given by 


2 | 
- Uy 
Cp =]- ta (25 ) 
The Ch change due to the increasing 61 is called the buoyancy correction. For the 


four stations where the boundary-layer profile was measured, the buoyancy- 
induced Cp was computed and is presented in Figure 29 as solid square symbols. 


The two curves shown provide reasonable agreement and illustrate why the 
pressure gradient was negative. | 


| | c-iiF-“buoyanc 
C | 3 : — icorrection 





Figure 29. Effect of buoyancy on pressure distribution on the ceiling of the 24” 
long test section. 
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b. Adverse Pressure Gradient 


The pressure distributions measured in the original test section and 
in the test section with the diffuser are presented in Figure 30. When either test 
section diffuser (7° or 14°) was installed the reference velocity (or dynamic 
pressure) was measured by a pitot-static probe located at station 1". The surface 
pressure data from both test section lengths showed a difference between the test 
section ceiling flat surface and the test section floor with a deflected surface. On 
the test section floor at station 6" (or 102"), the surface had an obtuse angle where 
there was an abrupt change of flow direction. As discussed in Milne-Thompson 
[Ref. 87] for potential flow, the velocity is locally infinite at an angle apex. The 
effects of a sharp corner in viscous flow also show an increased, but finite, 
velocity. At the angle apex, measured pressure coefficients were of the order of 
minus one and increased rapidly a small distance from the corner to a few tenths. 
These flows are discussed in greater detail by Mason [Ref. 88]. The data from the 
present investigation confirms the trends discussed above. The floor pressure 
distribution had a discontinuity near the angle apex at station 6" (or 102") where 
the local pressure coefficient peaked as low as -0.19 at x = 7.4" for the 24" test 
section and at about -0.1 at x = 5" for the 120" test section. These pressure 
coefficients correspond to 10% or 5% velocity increases, respectively. 


The inviscid pressure distribution for both test section diffusers was 
computed using the low-order panel method PMARC [Ref. 89] and the results for 
the 7° diffuser are presented in Figure 31. Without viscosity the pressure peak 
near the bottom apex angle is seen to peak at about -0.10 for the 7° diffuser which 
is similar to the experimental data and at about -1.40 for the 14° diffuser which is 
much lower than the experimental data in Figure 31. This indicates that while the 
flow for the 7° diffuser was attached, the flow for the 14° diffuser may have 
separated near the angle apex and as a result generate a reduced pressure peak. 
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Figure 30. Measured longitudinal pressure gradient. 
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(a) SJSU 24" test section. 
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(b) SJSU 120" test section. 
Figure 30. Concluded 
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Figure 31. Test-section pressure distribution computed by the inviscid panel- 
method PMARC [Ref. 90] for the 24" test-section with the 7° test-section diffuser. 


Since the velocity was subsonic, the streamlines above the angle 
apex should quickly smooth out, and the disturbance die down exponentially with 
increasing distance [Ref. 88]. The ceiling pressure-distributions for all six test 
section configurations show the expected and desired smooth profile in the vicinity 
of the cavity location(s). The 14° diffuser data still has a slope change near station 
6" (or 102") and shows the effect on the ceiling pressure-distribution of the change 
in pressure gradient slope at station 6" (or 102"). The comparison between the two 
test sections shows that the thicker boundary layer (see Figure 30) in the 120" test 
section was more effective in smoothing the Cp distribution on the ceiling 
pressure-distribution. These data also show that there was no significant flow 
separation apparent in these experimentally measured profiles. 


The experimental pressure distributions in the 0°, 7°, and 14° 


diffuser constant-angle portions of both the 24" and 120" long test sections are 
presented in Figure 32. These data show a linear variation of Cp with x location. 


Repeatibility is shown by the data for the 14° diffuser. 
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(a) 24" long test section. 
Figure 32. Effect of diffuser angle on the test section ceiling pressure distribution. 
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(b) 120" long test section. 
Figure 32. Concluded. 
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The resultant pressure gradients are tabulated below and plotted in Figure 33. 


Test Section ' Diffuser Angle, height at test 
Length, in. section diffuser 
entrance(x=6"), in. 


-0.00064 . 
0.01103 





0.06 


— {]— 24" test section 


a= 120" test section 
0.04 | 


d(Cp)/dx 


0.02 





diffuser angle 


Figure 33. Test section pressure gradient as a function of the diffuser angle. 
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Also shown in Figure 33 is the pressure coefficient for inviscid flow. 
This curve was computed using equation 24 with U; = U,,, the velocity at the test 


_ section entrance, and equation 25. The effect of boundary layer growth is to 


increase local velocity above that computed inviscidly, as described above for 
Figure 29. The slopes shown in Figure 33 are listed below: 


(zfs 
dx | 


These data show a 41% and a 60% reduction for viscous effects on the pressure 
gradient per degree of diffuser angle in the 24" and 120" long test sections, 
respectively. 










120" long test section 





Boundary-layer profiles were measured on the ceiling of the test 
section at four longitudinal stations in each test section: x = 5.5", 9", 12", and 16" 


in the 24" test section; or x = 101.5", 105", 108", and 112" in the 120" test section. 


These boundary-layer profiles are presented in Appendix D in the following 
figures: 
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Configuration Test Section 
length, in 


1400 
8400 


These profiles were used to compute the boundary-layer thickness 46, displacement 
thickness 5,, momentum thickness 55, pressure-gradient parameter 6 , the ratio 





of Clauser defect thickness to boundary-layer thickness A/od, and shape factor G 
which are presented in the following figures: 


Parameter 





Figure 


The data from the 24" test section are presented in part (a) and from the 120" test 
section are presented in part (b) of each of the above figures. 


The boundary-layer thickness, boundary-layer displacement 
thickness, and the boundary-layer momentum thickness data for the basic test 
sections (configurations 0000 and 8000) without pressure gradients were 
compared in Figure 23 with results from application of the 1/7 power law. While 
the 24" test section data were found to agree with the 1/7 power law, the higher 
Reynolds number 120" test section data tended to agree with a lower order 1/11 
power law. For the 24" test section the addition of the 7° and 14° diffusers 
reduced the values of 5, 5;, and 65 at corresponding stations in the divergent 


section (x = 9", 12", and 16"). In contrast, for the 120" test section the addition of 
the 7° and 14° test-section diffusers increased the values of 6, 0;, and 65 at the 


corresponding stations in the divergent section (x = 105", 108", and 112"). 


The Clauser equilibrium parameter B is often used to characterize 
the magnitude of a pressure gradient. The dominant term in this parameter is the 
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pressure gradient, dp/dx. This term is multiplied by the ratio of displacement 
thickness to skin friction. In Clauser's paper [Ref. 8], constant values of B were 
used to characterize equilibrium boundary layers. The results in Figure 37 show 
that the present investigation achieved equilibrium for configurations 0000, 8000, 
and 0700. The Clauser equilibrium parameter for the 24" test section presented in 
Figure 37(a) shows a variation similar to the actual pressures shown in Figure 
32(a). The Clauser equilibrium parameter for the 120" test section presented in 
Figure 37(b) also shows a variation similar to the actual pressures shown in Figure 
32(b) except for the 7° diffuser. Inspection of the values that make up the Clauser 
equilibrium parameter show that the skin friction coefficients and displacement 
thicknesses were nearly constant, which lead to a small variation of B. 


Another measure of the pressure gradient is the Clauser shape factor 
[Ref. 8], which is the ratio of the second to the first moment about the axis (U- 
u)/u* = 0. Another useful parameter is the Clauser defect thickness A. The 
relation of G to the universal parameter A/S was originally shown by Clauser 
[Ref. 8]. Equivalent information from the present investigation are shown in 
Figure 38. These data indicate that the test section diffusers induced a relatively 
modest adverse pressure gradient. 


The conventional shape factor H_ is often used as an indication of 


flow separation. Schlichting [Ref. 4, page 630] indicates that separation occurs | 


when H = 1.8 to 2.4. The data in Figure 39 show that all of the test sections gave 
Shape factors which were well below 1.8. It is interesting to observe that the no- 
diffuser configuration for the 24" test section had shape factors greater than either 
of the two diffuser configurations. In contrast, the opposite relation was measured 
in the 120" test section. All of these data show that the goal of obtaining adverse 
pressure gradients without flow separation was achieved in the present 
investigation. | 


The log-law relation (equation 18) for the data from the SJSU low 
speed tunnel was presented in Figure 25 and showed the expected agreement. The 
corresponding data for the 7° and 14° test section diffusers are presented in Figure 
40. The data for the 24" test section (Figures 40(a) and (b)) show reasonable 
agreement with the law-of-the wall equation using the generally-accepted 
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empirical constants. In contrast, the 120" test-section data fell below the law-of- 
the wall equation. This is consistent with the fact that the boundary layer 
thickness 5 was less than that given by the 1/7 power law. It appears that this 
boundary layer thickness is characteristic of a higher Reynolds number flow 
consistent with a 1/9 to as low as a 1/11 power-law correlation depending on the x 
station of the boundary-layer survey. 


Clauser [Ref. 8] used a universal plot of equilibrium turbulent 
velocity profiles, (U-u)/u* versus y/A, to display the effect of boundary-layer run 
length on the boundary-layer shape. The increasing x distance was shown in 
Figure 37(a) to correspond to the greatest increase in the pressure-gradient 
parameter $B achieved in this investigation. The resultant effect of the pressure 
gradient on the velocity profile is shown in Figure 41(b) using data from the 24° 
test section with the 14° diffuser. The changing velocity profile from station 5.5" 
to station 16" due to increasing { is consistent with the changes shown by 


Clauser. 


Since these profiles seemed to depend more on x location than 
specific diffuser geometry, the data from Figure 41 was replotted at a given x 
location for the 24" test section in Figure 42 and for the 120" test section in Figure 
43. The data at station 5.5" (Figure 42(a)) showed an anomaly in that an 
apparently large effect of pressure gradient was present. For this station the 7° 
data was apparently invalid and the 14° data showed a large profile change. The 
- other three stations from the 24" test section showed only small profile changes 
consistent with the data correlation shown in Figure 38. For the 120" test section 
the expected larger profile changes due to the larger pressure gradient parameter 
variation are seen. At x = 101.5" there was no pressure gradient, so the boundary 
layer profiles are seen to be in equilibrium. That is, all three surveys are the same. 
At the next three stations there was an increasing gradient in the pressure variation 
with increasing x which is consistent with increasing pressure gradient parameter, 


B. 
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(a) 24" test section. 
Figure 34. Effect of test-section diffuser angle on the variation of boundary-layer 
thickness as a function of longitudinal location. 
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(b) 120" test section. 
Figure 34. Concluded. 
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(a) 24" test section. 
Figure 35. Effect of test-section diffuser angle on the variation of boundary-layer 
displacement thickness as a function of longitudinal location. 
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(b) 120" test section. 
Figure 35. Concluded. 
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(a) 24" test section. 
Figure 36. Effect of test-section diffuser angle on the variation of ee layer 
momentum thickness as a function of longitudinal location. 
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(b) 120" test section. 
Figure 36. Concluded. 
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(a) 24" test section. 
Figure 37. Effect of test-section diffuser angle on the variation of freestream 
-pressure-gradient parameter as a function of longitudinal location. 
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(b) 120" test section. 
Figure 37. Concluded. 
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(a) 24" test section. | 
Figure 38. Effect of test-section diffuser angle on the ratio of the Clauser defect 


thickness to the boundary-layer thickness as a function of the Clauser shape factor. 
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(b) 120" test section. 
Figure 38. Concluded. 
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| (a) 24" test section. 
Figure 39. Effect of test-section diffuser angle on the variation of shape factor as a 
function of longitudinal location. 
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(b) 120" test section. 
Figure 39. Concluded. 
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(a) 24" test section with 7° diffuser. 
Figure 40. Effect of test-section diffuser angle on the variation of boundary-layer 
thickness as a function of longitudinal location. 
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(b) 24" test section with 14° diffuser. 
Figure 40. Continued. 
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(c) 120" test section with 7° diffuser. 
Figure 40. Continued. 
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(d) 120" test section with 14° diffuser. 
Figure 40. Concluded. 
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(a) 24" test section with 7° diffuser. 
Figure 41. Effect of test-section survey longitudinal position on the universal 
velocity distribution. 
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(b) 24” test section with 14° diffuser. 
Figure 41. Continued. 
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(c) 120" test-section with 7° diffuser. 
Figure 41. Continued. 


(U-u)/u* 





(d) 120" test-section with 14° diffuser. 
Figure 41. Concluded. 
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(a) Station 5.5" in the 24" test section. 
Figure 42. Effect of test-section diffuser angle on the universal velocity profiles. 
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— (b) Station 9" in the 24" test section. 
Figure 42. Continued. 
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(c) Station 12" in the 24" test section. 
Figure 42. Continued. 
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(d) Station 16" in the 24" test section. 
Figure 42. Concluded. 
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| (a) Station 101.5" in the 120" test section. 
Figure 43. Effect of test-section diffuser angle on the universal velocity profiles. 
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(b) Station 105" in the 120" test section. 
Figure 43. Continued. 
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(c) Station 108" in the 120" test section. 
Figure 43. Continued. 
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(d) Station 112" in the 120" test section. 
Figure 43. Concluded. 
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2. Open Cavity Configurations 


The purpose of the present investigation was to experimentally and 
computationally evaluate the effect of 1, 2, or 4 two-dimensional spanwise cavities 
with square cross-sections on the boundary-layer characteristics in both a constant 
pressure flow and two adverse pressure-gradient flows. The experimental results 
presented in the present section are based on the boundary-layer profiles presented 
in Appendix D and the surface/cavity pressure distributions presented in Appendix 
C. These data will be compared and analyzed in the following two subsections. 


a. Effects on the Boundary-Layer Characteristics 


The effects of the cavity geometry on the boundary-layer 
characteristics are analyzed in this subsection to identify how they are influenced 
by the flowfield. The boundary layer profiles measured on the tunnel ceiling in 
the vicinity of the cavities were used to compute the boundary layer thickness 6, 
displacement thickness 6, (equation 4), momentum thickness 05, shape factor H, 
and skin-friction coefficient cr which are presented in the following figures: 





Figure 


The data from the 120" test section are presented in parts (a), (b), and (c) (test 
section diffuser angles of 0°, 14°, and 7°, respectively) and from the 24" test 
section are presented in parts (d), (e), and (f) (test section diffuser angles of 0°, 
(14°, and 7°, respectively) of each of the above figures. 


The longer 120" test section provided a fully-developed turbulent 
boundary layer with only small increases in the x = 100" to 112" measurement 


region where the Reynolds number was large (Re > 5x10®). The boundary-layer 
thicknesses (0, 61, and 05) for the 120" test section with either the 0° or the 7° test. 
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section diffusers did not change significantly with variation in the cavity 
configuration (parts (a) and (c) of Figures 44 through 46). Depending on the 
selected comparison, these data were within + 2% to + 12% for a given plot. Only 
the four-cavity configuration with the 14° diffuser shows a noticeable increase (up 
to 50%) in these thicknesses (part (b) of Figures 44 through 46). 


| The 24" test section had a much thinner boundary layer which is 
initially barely turbulent since the Reynolds number ranged from 350,000 to about 
1,100,000. The boundary layer profiles presented in Appendix D tend to agree 
with the 1/7 power law which is typical of a two-dimensional turbulent boundary 
layer. The boundary layer thicknesses (6, 01, and 6) for the 24" test section with 
either the 0° or 7° test section diffuser are typically within + 10%. In the 14° test 
section diffuser, while there was more data scatter, there tended to be an increase 
in thickness with the increasing diffuser angle. The shape factor H (Figure 47) 
and skin friction coefficient c¢ (Figure 48) do not show any noticeable effect of 


the presence of the cavities. 


The boundary layer log-law profiles were compared (Figures 49 to 
51) both upstream and downstream of the cavity (or cavities) in an attempt to 
identify their effect on the flow. It is shown that all of the data obtained without a 
diffuser tend to agree with the log law (equation 18) where « = 0.41 and Cj = 5.0. 


The data with the 7° and 14° test-section diffusers asymptote below the log-law 
defined by the above constants. It was found, as shown in Figure 52, that if Cj is 
changed to 3.5 for the 7° diffuser and to 3.0 for the 14° diffuser that these data for 
the test-section diffuser geometries agree with the log-law equation. 


The effect of variation of the cavity configuration (0, 1, 2, or 4 
cavities) on the boundary-layer profiles at a given x station is presented in Figures 
53, 54, and 55. The test-section entrance profiles (Figure 53) demonstrate 
- boundary layer similarity for the basic test section and for the 7° test section 

diffuser. These profiles are in agreement with the 1/7 power law profile. There is 
a significant boundary layer profile variation from the 1/7 power law for the 14° 
test-section diffuser. The data from Figure 54, which are for the downstream x 
station, are presented in Figure 55 using the universal velocity profile. This data 
format demonstrates a modest variation in the boundary-layer profiles which are 
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typical of an adverse pressure gradient. The single and double cavities show only 
a small deviation from the profile measured with no cavity. The four-cavity 
configuration profile in most cases represents a more adverse boundary-layer 
profile than any other configuration. 
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b. Effects on the Surface Pressure Distribution 


The effect of the test-section diffuser on the surface-pressure 
distribution was discussed on page 82. The focus of the present section is the 
evaluation of the effect of cavities on the pressure distribution on the wind-tunnel 
ceiling. Several configurations were measured several times to determine the data 
repeatibility. An example is presented in Figure 56 where the ceiling pressure 
distributions from five runs of configuration 8431c (120" test section with a 2" 
deep cavity covered by a door located between x = 102" and 104") are plotted 
along with their average. Immediately downstream of the cavity the pressure 
coefficient is seen to be depressed. This was probably caused by a misalignment 
of the door with the ceiling. Most of the pressure coefficients are within + 0.02 of 
the average pressure coefficient curve shown on Figure 56. 


The effect of the cavities in the 24" test section is presented in 
Figures 57 and 58. The basic test section shows a pressure coefficient magnitude 
of approximately -0.02; ahead of the cavities there is a similar pressure coefficient. 
Downstream of the cavity the pressure coefficients drop to about -0.06. With both 
test section diffusers (Figures 59 through 61) in the 24” test section the general 
data trend is unaffected by the presence of the cavities. 


The effect of the cavities in the 120" test section is presented in 
Figures 62 through 66. Generally there is seen to be little effect of the cavities on 
the pressure distributions. The one exception is the effect of the 2" cavity 
downstream of that cavity where the pressure coefficient is reduced to about -0.04. 
' With both test section diffusers there is more scatter in the data without a 
consistent trend. 


In conclusion, all of the surface pressure data suggest that the 
presence of the cavities has no consistent effect. It may be concluded that cavities 
may be deployed without significantly modifying the resultant pressure 
distribution. 
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(a) 120" test section. 
Figure 44. Effect of the cavity configuration on the boundary-layer thickness. 
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(b) 120" test section with the 14° test-section diffuser. 
Figure 44. Continued. 
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(c) 120" test section with the 7° test-section diffuser. 
Figure 44. Continued. 
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(d) 24" test section 
Figure 44. Continued. 
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(e) 24" test section with the 14° test-section diffuser. 
Figure 44. Continued. 





(f) 24" test section with the 7° test-section diffuser. 
Figure 44. Concluded. 
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| (a) 120" test section. | | 
Figure 45. Effect of the cavity configuration on the displacement thickness. » 
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(b) 120" test section with the 14° test-section diffuser. 
Figure 13. Continued. 
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(c) 120" test section with the 7° test-section diffuser. 
Figure 45. Continued. | 
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(d) 24" test section 
Figure 45. Continued. 
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(e) 24" test section with the 14° test-section diffuser. 
Figure 45. Continued. 
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(f) 24" test section with the 7° test section diffuser. 
Figure 45. Concluded. 
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(a) 120" test section. 
Figure 46. Effect of the cavity configuration on the momentum thickness. 
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(b) 120" test section with the 14° test-section diffuser. 
Figure 46. Continued. 
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(c) 120" test section with the 7° test-section diffuser. 
Figure 46. Continued. 


no. cavities, D, in. 
—O-—— 0, = O=- 2, 1 
wegen 1, 2 wmmp\ume 4, | 





X, in. 


(d) 24" test section 
Figure 46. Continued. 
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(e) 24" test section with the 14° test-section diffuser. 
Figure 46. Continued. | 
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(f) 24" test section with the 7° test-section diffuser. 
Figure 46. Concluded. 
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(a) 120" test section. 
Figure 47. Effect of the cavity configuration on the shape factor. 
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(b) 120" test section with the 14° test-section diffuser. 
Figure 47. Continued. 
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(c) 120" test section with the 7° test-section diffuser. 
Figure 47. Continued. 
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(d) 24" test section 
Figure 47. Continued. 
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(e) 24" test section with the 14° test-section diffuser. 
Figure 47. Continued. | 
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(f) 24” test section with the 7° test-section diffuser. 
Figure 47. Concluded. 
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(a) 120" test section. 


Effect of the cavity configuration on the skin-friction coefficient. 


0.006 ——— 
no. cavities, ‘D, in. 


0.004 


0.002 





0 
100 104 108 112 
X, in. 


(b) 120" test section with the 14° test-section diffuser. 
Figure 48. Continued. 
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(c) 120" test section with the 7° test-section diffuser. 
Figure 48. Continued. 
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(d) 24" test section 
Figure 48. Continued. 
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(e) 24" test section with the 14° test-section diffuser. ' 
Figure 48. Continued. 
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(f) 24" test section with the 7° test-section diffuser. 
Figure 48. Concluded. 
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(a) 120" test section at x ~ 101". 
Figure 49. Effect of the cavity geometry on the boundary-layer log-law profile. 
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(b) 120" test section at x ~ 101" with the 14° test-section diffuser. 
Figure 49. Continued. 
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(c) 120" test section at x ~ 101" with the 7° test-section diffuser. 
Figure 49. Continued. 
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(d) 24" test section at x = 5". 
Figure 49. Continued. 
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(e) 24" test section at x ~ 5" with the 14° test-section diffuser. 
Figure 49. Concluded. | 
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(a) 120" test section at x =~ 108". 
Figure 50. Effect of the cavity geometry on the boundary-layer log-law profile. 
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(b) 120" test section at x ~ 108" with the 14° test-section diffuser. 
Figure 50. Continued. 
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(c) 120" test section at x ~ 108" with the 7° test-section diffuser. 
Figure 50. Continued. 
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(d) 24" test section at x = 12”. 
Figure 50. Continued. 
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(e) 24" test section at x ~ 12" with the 14° test-section diffuser. 
Figure 50. Continued. 
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(f) 24" test section at x = 12" with the 7° test-section diffuser. 
Figure 50. Concluded. 
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(a) 120" test section at x= 112". | 
Figure 51. Effect.of the cavity geometry on the boundary-layer log-law profile. 
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(b) 120" test section at x = 112" with the 14° test-section diffuser. 
| Figure 51. Continued. 
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(c) 120" test section at x = 112" with the 7° test-section diffuser. 
| Figure 51. Continued. | 
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(d) 24" test section at x = 16". 
Figure 51. Continued. — 
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—(e) 24" test section at x = 16" with the 14° test-section diffuser. 
Figure 51. Continued. 
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(f) 24" test section at x = 16" with the 7° test-section diffuser. 
Figure 51. Concluded. 
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(a) 7° test-section diffuser data compared with log-law equation using C] = 3.5. 


Figure 52. Comparison of 120" test-section experimental data at x = 112" with the 


log-law equation with C; changed from 5.0. 
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(b) 14° test-section diffuser data compared with log-law equation using C, = 3.0. 
Figure 52. Concluded 
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(a) 120" test section at x = 101”. 
Figure 53. Effect of the cavity configuration on the boundary-layer profile. 
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(b) 120" test section at x = 101" with the 14° test-section diffuser. 
Figure 53. Continued. 
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(c) 120" test section at x = 101" with the 7° test-section diffuser. 
Figure 53. Continued. 
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(d) 24" test section at x = 5". 
Figure 53. Continued. 
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(e) 24" test section at x ~ 5" with the 14° test-section diffuser. 
_ Figure 53. Concluded. 
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(a) 120" test section at x = 112”. 
Figure 54. Effect of the cavity configuration on the boundary-layer profile. 
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(b) 120" test section at x = 112" with the 14° test-section diffuser. 
Figure 54. Continued. 
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©) 120" test section at x = 112" with the 7° test-section diffuser. 
Figure 54. Continued. 
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(d) 24" test section at x = 16". 
Figure 54. Continued. 
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(ec) 24" test section at x = 16" with the 14° test-section diffuser. 
Figure 54. Continued. 
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(f) 24" test section at x = 16" with the 7° test-section diffuser. 
Figure 54. Concluded. 
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(a) 120" test section at x = 112". 
Figure 55. Effect of the cavity configuration on the universal velocity profile. 
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(b) 120" test section at x = 112" with the 14° test-section diffuser. 
Figure 55. Continued. 
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(c) 120" test section at x = 112 with the 7° test-section diffuser. 
Figure 55. Continued. 
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(d) 24" test section at x = 16". 
Figure 55. Continued... 
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(e) 24" test section at x = 16" with the 14° test-section diffuser. 
Figure 55. Continued. 
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(f) 24" test section at x = 16" with the 7° test-section diffuser. 
Figure 55. Concluded. 
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Figure 56. An example of pressure-data repeatability using the tunnel 
ceiling measurements for configuration 8431c (120" test section and D = 2" with 
the door closed). 





Figure 57. Effect of the single cavity size on the ceiling pressure distribution (24" 
_ test section with D = 1" or 2"). 


142 








Figure 58. Effect of the multiple cavities on the ceiling pressure distribution (24" 
| test section with D = 1"). 





Figure 59. Effect of the single cavity size on the ceiling pressure distribution (24" 
test section, 14° test-section diffuser with D = 1" or 2"). 
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Figure 60. Effect of the multiple cavities on the ceiling pressure distribution (24° 
test section, 14° test-section diffuser with D = 1"). 





Figure 61. Effect of the multiple cavities on the ceiling pressure distribution (24° 
test section, 7° test-section diffuser with D = 1"). 
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Figure 62. Effect of the single cavity size on the ceiling pressure distribution 
(120" test section with D = 1" or 2"). 





Figure 63. Effect of the multiple cavities on the ceiling pressure distribution (120" 
test section with D = 1"). 
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Figure 64. Effect of the single cavity size on the ceiling pressure distribution 
(120" test section, 14° test-section diffuser with D = 1" or 2"). 
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Figure 65. Effect of the multiple cavities on the ceiling pressure distribution (120° | 
test section, 14° test-section diffuser with D = 1"). 
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Figure 66. Effect of the multiple cavities on the ceiling pressure distribution (120" 
test section, 7° test-section diffuser with D = 1"). 


3. Cavity With a Door 


As described in the introduction, the effect of cavities are being studied as a 
means for increasing airfoil camber to provide a new high lift device. One aspect 
of this application is the effect of surface non-uniformities (which could be caused 
by cavity doors) on the boundary layer and on the surface pressure distribution. 
Tests were conducted which compared a plain, smooth ceiling with a ceiling 
which had a two-inch square cavity which was closed by a flush door. The door 
was mounted by a pivot at its upstream end. When the wind tunnel was not 
operating, the closed door was adjusted to be flush with the adjacent tunnel 
ceiling. When the tunnel was turned on and operated at a freestream velocity of 
120 feet per second, it was noted that the door deflected into the flow 
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approximately 1/16 to 1/8 inch at its downstream edge. The larger deflections 
occurred when one of the test-section diffusers was present ( @ = 7° or 14°). 


For these tunnel geometries, the measured boundary-layer thicknesses are 
presented in Figure 67 and the pressure distributions are presented in Figure 68. In 
each figure, part (a) presents data from the 24" test section and part (b) has data 
from the 120" test section. The data for the plain ceiling are presented by the solid 
symbols and the data for the cavity with a door are presented by open symbols. 
The boundary-layer thicknesses for the 24" test section cavity with a diffuser and a 
cavity with a closed door (Figure 67(a)) are thicker than those for the plain ceiling. 
The boundary layer thickness for the closed cavity with the 7° test-section diffuser 
increased much more rapidly than any of the other configurations. This difference 
is probably due to the increased projection of the cavity door trailing edge into the 
freestream, which was a significant fraction of the 6 of 0.25 to 0.45 inch in this 
test section. The data for the 120" test section (Figure 67(b)) show a more 
noticeable effect of the cavity door projection even though the door projection was 
a small fraction of the measured 6, which was greater than an inch. This indicates 
that poor fit of a cavity door could increase skin friction drag. It is seen in Figure 
68 that there were differences in the longitudinal pressure distribution for both test 
sections. This indicates that the surface non-uniformity had an effect on the 
inviscid freestream which determines the surface pressure distribution. 


A second set of tests was conducted which compared the flow in the 
- presence of a two-inch square open cavity with a similar two inch cavity ( the 
closed cavity used in Figures 67 and 68) where a door was opened into the cavity 
to form a front cavity wall of an open-door cavity. The open cavity without a door 
had a smooth surface upstream of the cavity. In contrast, the open cavity door had 
a small gap of about 1/16 inch from the upstream end of the ceiling surface. These 
results are presented in Figure 69 (boundary layer thicknesses) and Figure 70 
(longitudinal pressure distribution). While these results appear to be similar to 
those presented for the closed cavity in Figures 67 and 68, there are notable 
differences. The boundary-layer thicknesses for all of the open-door cavity 
configurations showed an increase with longitudinal distance similar to the — 
increase of the closed cavity which had the largest increase, the 7° test-section 
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(a) 24" test section 
Figure 67. Effects of closed cavity door and test-section diffusers on the variation 
of boundary-layer thickness as a function of longitudinal location. 
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(b) 120" test section 
Figure 67. Concluded 
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(a) 24" test section 
Figure 68. Effects of closed cavity door and test-section diffusers on the variation 
of ceiling pressure distribution. 
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(b) 120" test section 
Figure 68. Concluded 
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(a) 24" test section 
Figure 69. Effect of open cavity door and test-section diffuser angle on the ceiling 
boundary-layer thickness as a function of longitudinal location. 
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(b) 120" test section 
Figure 69. Concluded 
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(a) 24" test section 
Figure 70. Effect of open cavity configuration and test-section diffuser angle on 
the ceiling longitudinal pressure distribution. 
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(b) 120" test section 
Figure 70. Concluded 
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configuration. The reasons for this result have not been identified. In contrast, it 


is seen that with the open cavities (Figure 70) there were smaller differences in the 


longitudinal pressure distribution for either test section length than that shown by 


the closed cavity and plain ceiling. These results suggest that the open cavities 


produce less change in the surface pressure distributions than do surface non- 
uniformities in an otherwise plain ceiling. 
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IV. INCOMPRESSIBLE NAVIER-STOKES COMPUTATIONS 


There are some simple flows related to the present investigation which may 
be represented by empirical correlations. One example is a two-dimensional, 


turbulent boundary-layer profiles using the 1/7 power law (u/U,, =(y/8)\" 7) 


More complete flow solutions require more complete sets of equations. There is a 
range of equation sets of increasing complexity available for the solution of fluid 
flows. All of these equation sets may be obtained as simplifications of the Navier- 
Stokes equations. The particular equations selected for a problem of interest 
depends on the characteristics of that flow. The simplest of these equations is the 
linear potential Laplace equation which assumes inviscid, incompressible, 
irrotational flow. 


V70=0 (26) 


When there are viscous effects more complex equation sets are needed. Two- 
dimensional, viscous flows in a conventional, subsonic wind tunnel test section 
may usually be accurately computed using the boundary-layer equations 
(Equations 11 through 13). These equations describe the attached boundary layer 
near the walls which develops adjacent to the inviscid outer flow, which is 
representative of the bulk of the freestream flow. 


The boundary-layer concept of Prandtl assumes that the viscous layer is 
small (thin boundary layer) relative to the streamwise flow distance (6/L <<1). 
As a consequence, the Navier-Stokes equations may be simplified to obtain the 
boundary layer approximation for steady, two-dimensional, incompressible, 
constant fluid-property flow [Ref. 90]. The present experimental investigation 
adds an adverse pressure gradient and one or more cavities to produce a flow 
which is more complicated than those which can be resolved by the boundary- 
layer equations. Especially difficult is the resolution of the recirculating flow in a 
cavity and its interaction with the freestream in an adverse pressure gradient. For 
example, in the 120-inch test section in the present investigation, the cavity flow 
had a viscous layer whose thickness was nearly the same as the cavity depth. In 
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contrast, for the 24 inch test section the viscous boundary layer thickness was 
much less than the cavity depth. 


Based on the present literature survey, neither experimental nor 
computational data have been obtained for these flows prior to the present 
investigation. The present experiment measured the flow in the vicinity of several 
cavity configurations with a freestream Mach number of the order of 0.10. For 
flows such as these, the incompressible Navier-Stokes equations are a more 
suitable representation than the boundary layer equations. The incompressible 
Navier-Stokes equations are derived from the Reynolds-averaged Navier-Stokes 
equations by assuming an incompressible flow (i.e., M = 0 and a= 9). 


Since the temperature was essentially constant in the present experimental 
investigation the energy equation is not required. The two-dimensional, 
incompressible Navier-Stokes equations in Cartesian coordinates become the 
following: 


continuity equation, ——+—=0 (27) 


X momentum equation OF gt du, au (28) 
: "ot ax ody pox 2 oy” } 


) 
y-momentum equation, —+u—+Vv—=--— 


A numerical solution of these equations may be obtained for example using either 
a finite-element or finite-difference scheme, a suitable grid, and the definition of 
appropriate boundary conditions. 
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A. INCOMPRESSIBLE NAVIER-STOKES CODE (ins2d) 


The ins2d computer code [Ref. 91 and 92] numerically solves the two- 
dimensional, incompressible Navier-Stokes equations in a time-accurate manner. 
In a pseudocompressibility formulation, a time derivative of pressure is added to 
the continuity equation, 


ou ov 
22 +{ Ss +e 0 (30) 


Note that B is the pseudocompressibility constant and t represents a pseudotime 
which is not related in any way to physical time. Equation 30 is then combined 
with the momentum equations to obtain the following equations in two- 
dimensional Cartesian coordinates. 


) e) e) 
-——D E-E, )+—(F-F,})=0 31 
iad Plies |e (31) 
where 
Dp Bu Bv 
D=lu E= u-+p Fe uv 
uy V +p 
0 0 
Ey =] Tx FY = Tyx 
txy Tyy 
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The convective terms are upwind-differenced using a Roe [Ref. 93] flux- 
difference split approach that has uniformly high accuracy throughout the interior 
grid points. This scheme was derived as an approximate Riemann solver for the 
compressible flow equations. In this formulation the Reynolds stress has been 
approximated as a function of the strain-rate tensor, and thus represents a sum of 
the kinematic viscosity and the turbulent eddy viscosity. Constant kinematic 
viscosity is assumed for simplicity. The viscous fluxes are differenced using 
second-order accurate central differences. 


The system of equations (equations 31) is solved numerically using an 
implicit line-relaxation scheme. Application of a first-order backward Euler 
formula to equations 31 yields the delta-form equation 


I oR ‘ n+i pyn\__pa 
ee(Z) (> D )= R (32) 


Here the superscript n is the pseudotime iteration count and R is the residual 
vector. Unsteady problems are solved with the use of subiterations in pseudotime 
at each physical time step. 


At the inflow boundary there is one characteristic wave traveling out of the 
computational domain the since fluid is traveling into the domain. When the 
inflow velocity profile is not known, there is a vector of variables which is held 
constant and defined as 
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uv 
0 0 (33) 
0 1 


For the present calculation the vertical inflow velocity component v was set to 
zero and the horizontal inflow velocity u was set to unity at all of the inflow 
plane grid points except at the floor and ceiling grid points where both were zero. 


At the outflow boundary there are two characteristic waves traveling out of 
the computational domain since fluid is also leaving the domain. This 
computation used a specified static pressure at the outflow plane which was 
defined as 


p 1 O Q 
o=! 0 o lo 0 0 (34) 
aD 
0 0 0 


The tangential flow boundary conditions on the tunnel floor and ceiling were 
obtained by setting both the u and v velocity components to zero. 


B. GRID GENERATION 


The grid was generated using several computers and codes. Several 
FORTRAN codes were written to generate both the freestream flow grid between 
the wind tunnel ceiling and the floor (sjtgg.f and rtgg.f) and the cavity grid (cgg.f 
and cgga.f). Then these grids were coupled by the PEGSUS code [Ref. 94] which 
uses an overset scheme. Typically these codes were run interactively. For the 
basic wind tunnel flow computations the grid from code sjtgg.f was used directly 
in the ins2d code. These FORTRAN programs were solved using a Silicon 
Graphics (SGI) Challenger L workstation. This computer has 4 processors, a 
MIPS R4400 microprocessors with a clock speed of 150 MHz, 256 MB random- 
“access memory and 2GB random-access disk storage. 
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For a boundary layer flow, the largest velocity gradient occurs adjacent to 
the walls and normal to the flow direction. For the present ins2d Navier-Stokes 
flow solver, it is recommended that the finest spacing near the wall lead to a 
minimum non-dimensional inner layer distance normal to the wall yt of 
approximately 1. The yt is reduced by reducing the grid spacing near the tunnel 
wall. In Anderson, Tannehill, and Pletcher [Ref. 90] several simple 
transformations are presented which can be used to cluster the grid in regions of 
large gradients such as boundary layers. These transformations are part of a 
general family proposed by Roberts [see Ref. 90, pp. 247-250]. One 
transformation was developed to refine the mesh near the walls of a two- 
dimensional tunnel using equations 

X=X 
and 


y-a —O (35) 
aaenfie[+n/(-—} i 


y=h 


where h is tunnel height, o = 1/2 refines the grid equally near y = 0 and y =h, and 


h 
Coordinates X and ¥ represent a uniform Cartesian grid and coordinates x and y 


V7) 
p-(1-2 for 0<2 <I (36) 


represent a clusted grid. A computer program (rtgg.f) was written to generate the 
grid for the wind tunnel used in the experimental investigation. 


An example of Roberts’ transformation grid clustering is shown in Figure 
71(a) for the region between the tunnel floor and the tunnel centerline. The region 
from the tunnel centerline to the tunnel ceiling is the reflection of the clustering 
shown in Figure 71. The horizontal scale is a uniformly-spaced coordinate 1j 
non-dimensionized by the scale factor SF, and the vertical scale is clustered grid 
scale yj norm which is then converted from the 0 to 1 scale to the dimensions of 
the tunnel half height. The grid clustering is determined by the value of B which is 
shown to range from 2, for the least clustering, to 1.001, for the greatest clustering. 
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The largest two values, 2 and 1.2, provide only a limited clustering. In the next 
section results from the ins2d code will be presented for the three finest scale 
factors shown in Figure 71(a): 1.05, 1.01, and 1.001. These results will be used to 
evaluate the usefulness of the Roberts grid-transformation scheme for —" 
these boundary-layer flows. 


An alternative to Roberts’ transformation-clustering scheme was developed 
by the author using a tangential grid-spacing scheme. This transcendental function 
was used to transform a uniform spacing to a variable grid spacing. The fineness 
of the grid spacing was determined by a tangential grid-spacing scale factor SF. 
Initially the scheme divided a distance into a uniform spacing where 


5 max ~ Tmax ~ 1 a7 
N= Nex] (37) 

using 
nN, =tan'(-SF) and 7n,,,, =tan7!(SF) (38) 


The value of the scale factor is increased when a finer grid is desired at the outer 
grid edges. For the present investigation, SF was typically 2 in the streamwise 
direction and ranged from 16 to 72 normal to the freestream direction. 


As an illustration of the final steps only the y coordinate equations are used. 
Similar calculations are done in the x direction. For the next step, a normalized 
Spacing from -1 to 1 was computed as 


—| 
= tan 
yi = _tan Ti (39) 


—l 
tan “ Nmax 


Finally the physical coordinates were computed as 
= Vou + yi L 40 
yj org > ( | ) 


The fraction in the second term scales length from y; to a range from 0 to 1. 
Since Yorg is the minimum physical dimension and L is the dimension length, 
the yj values cover the entire length of this direction. Program sjtgg.f was used 
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to algebraically calculate a two-dimensional grid for the wind tunnel freestream 
including the tunnel ceiling and floor, with or without the test section diffuser. 
Program cgg.f was used to generate each cavity grid. Both calculations clustered 
the grid using a tangential spacing scheme. 


The normalized spacings obtained from this procedure is illustrated for 
several scale factors in Figure 71 (b). The increased clustering near the walls is 
readily apparent. At the highest scale factors the grid tends to be sparse near the 
center of the physical space. For the streamwise direction the flow tends to a zero 
gradient in the potential flow region away from the boundary layer. As a 
consequence the coarser spacing is adequate near the tunnel centerline. The grid 
size variation is well behaved for the present calculations because it changes by 
less than 1.2 between adjacent grid cells. This grid clustering scheme is evaluated 
in the next section to determine its suitability for the present investigation. 


0.75 


--O--1.05 
== 1.01 
---§--- 1.001 


y; norm 9.5 


0.25 





n/SE 
(a) Roberts grid transformation. 
Figure 71. Clustered grid spacing as a function of a uniform spacing for several 
scale factors. 
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(b) Tangential grid transformation. 
Figure 71. Concluded 


The tangential grid spacing is used for both directions in each cavity and for 
both tunnel directions. When the tunnel has one or more cavities, the tangential 
grid spacing in the streamwise direction must be modified to satisfy an adjacent 
grid spacing variation of 1.2 or less. This requirement was satisfied as part of the 
program cgga.f. First this program generates the grid for the specified cavity. 
Then the tunnel grid generated in sjtgg.f is read and modified. In the streamwise 
direction the grid spacing is reduced immediately upstream and downstream of the 
cavity. Generally the two grids, freestream and cavity, are compatibily spaced in 
the vertical direction and require no modification. If modification is needed to 
satisfy the 1.2 adjacent spacing criterion the program cgga.f accomplishes this 
change in the freestream grid. If more than one cavity is present then cgga.f 
continues on to add these cavities to complete the grid. 
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C. CFD RESULTS 


After the grid generation was completed, the CFD code computed the fluid 
flow solution and analysis of the results. After the grids were prepared, two 
additional input files were prepared which provided: (1) the boundary conditions 
and (2) the ins2d input data for the desired code options. The present input 
boundary conditions were (a) an inflow boundary with a constant velocity and 
characteristic relation for pressure, (b) an outflow boundary using characteristic 
relations for velocity and constant static pressure, (c) the ceiling with a no-slip 
wall with the wall normal vector pointing in the negative computational direction, 
and (d) the floor with a no-slip wall with the wall normal vector pointing in the 
positive computational direction. About 40 items of data make up the ins2d.in 
file. These data are defined in a user file described by the ins2d code and will not 
be discussed in detail in this paper. 


The ins2d code was run on a Cray C-90 supercomputer with 8 processors, 
256 Mw of random access memory (ram), 8 Gw of disk storage. The floating- 
point processor speed was 1 GFLOP which is about 20 faster than the Challenge L 
workstation used for grid generation. A typical case with a 201 by 321 grid used 
about 7.75 Mw of ram and required about 2.2 seconds/iteration or 345 x 10-/ 
seconds/iteration/grid point. The convergence criterion is based on the magnitude 
of the change in the residual term from the previous iteration. It was found that if 
a case converged it took 40 to 150 iterations. If a case took more iterations it did 
not converge. Some cases were run for as many as 20,000 iterations without 
convergence. 


For data analysis, the desired ins2d results were then selected from the 
solution file using a postprocessing program called plot3d on the Cray C-90. 
These results were then transferred as a group of ASCII data files using a file 
transfer protocol to an Apple Macintosh IIci personal computer. This computer 
used a Motorola 68030 microprocessor which had a clock speed of 40 MHz and an 
8 MB random-access memory and an 80 MB hard disk. In most cases these data 
were read by a plotting program called Kaleidagraph. 
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1. Wind Tunnel Flow Computation 


Several grid geometries were examined for the basic 24" test section 
(identified as C0000) and for the 120" test section (identified as C8000) used in 
the experimental investigation. For this phase of the computational investigation 
both the Robert's transformation and the tangential transformation were used to 
generate grids to represent the flow in the wind tunnel for both the 24-inch test 
section and the 120-inch test section. Also examined for both test section lengths 
were the number of grid points in the y direction, normal to the freestream. For 
the 24-inch test section, the effect of the longitudinal length of the x grid was 
varied. 


a. 24 Inch Test Section 


The grid was used to represent the 24 inch wind tunnel test section 
(C0000) from the floor to the ceiling and between the longitudinal stations of -16 
to 24 inches. A few preliminary cases were run for a shorter longitudinal grid 
which ran from 0 to 24 inches. It was found that this provided boundary layer 
thicknesses which were too thin in relation to the experimental data. The grid 
clustering for Roberts' transformation is determined by the value used for the 
stretching parameter $ (equation 36). For the 24" test section the boundary layer 
thickness 6 ranged from 0.22 to 0.45 inch depending on the longitudinal survey 
location. This lead to a recommended stretching parameter of about 1.01. In the 
ins2d code the resultant grid resulted in a minimum y+ value of 10.0 which is 
larger than the recommended value of about 1. | 


Using 201 grid points in the y direction, normal to the freestream, 
three values of B were used which resulted in a reduction in the minimum yt from 


14.0 to 1.8. An additional grid with 401 y grid points was also used to obtain a y+ 
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value near the desired value of 1. The following table outlines the Robert's 
transformation cases which were used to evaluate these grids. 


C0000 y grid points Robert's grid B 


a sore 


The ins2d code iterates until either the maximum number of iterations are 
completed or the computation converges. None of these cases converged within 
5000 iterations. This undesirable result indicates Lax's equilibrium theorem [Ref. 
90] is not being satisfied. This theorem states that given a properly-posed initial- 
value problem and a finite-difference approximation that satisfies the consistency 
condition, stability is the necessary and sufficient condition for convergence. 

The computed boundary-layer thicknesses (6) are presented in 









Figure 72 along with the experimental data and the 1/7 power law result. The 
latter two values are shown to be in good agreement. The computational result 
- from the coarsest grids produced thicknesses greater than experimental data. The 
computational result from the two finest grids where yt is near 1 show thicknesses 
which are less than the experimental data. These inadequate results reflect the lack 
of convergence in these computations and suggest that Roberts' transformation 
grids do not provide a satisfactory grid for the present problem. Even though the 
thicknesses do not show good agreement with the experimental data, the 
boundary-layer profiles at longitudinal locations 5.5 (Figure 73(a)) and 16 inches 
(Figure 73(b)) show agreement between the experimental and the computational 
results. This suggests that the profile is a less useful parameter for evaluating the 
adequacy of the computational results. 
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Figure 72. Boundary-layer thickness as a function of longitudinal tunnel location 
computed by the ins2d for 24" tunnel using Roberts’ transformation grid. 


1/7 power law : 
w= m= ins2d B=1.001 
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(a) x=5.3": 
Figure 73. Boundary-layer profile computed by the ins2d for 24" tunnel using 
Roberts’ transformation grid. 
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1/7 power law 2 
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(b) x = 16". 
Figure 73. Concluded. 


The tangential-grid scheme was also used for the 24 inch test section 
tunnel (C0000). Scale factors of 32 and 64 were selected to achieve yt values 
which approached 1. The number of y grid points was increased from 201 to 401 
to obtain a minimum value of 1.2. In contrast to the Robert's transformation all of 
these cases converged in from 48 to 76 iterations. These rapid solutions enabled 
the cases to be run in the debug queue on the Cray C-90 computer in less than 
three minutes. This rapid convergence provides confidence that these grids 
provide consistent, stable computational results. 






C0000 y grid points 


scale factor __ y+/ (iterations for convergence) __ 
tm 
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Examples of the computed boundary layer thicknesses using the tangential-grid 
scheme are presented in Figures 74(a) for 201 y grid-points and in Figure 74(b) 
for 401 y grid-points. Both figures present data for scale factors of 32 and 64. In 
both figures there is excellent agreement between the computational results, the 
experimental data, and the 1/7 power law at longitudinal locations at 5.5 inches 
and 12 inches. The reasons for the disagreement at the 16 inch station is 
unknown. It is possible there may be an influence of the diffuser downstream of 
the test section on the downstream end of the test section. This diffuser was not 
modelled in these calculations. Figure 75 presents comparisons of computational 
results, the 1/7 power law, and the experimental data for the boundary layer 
profiles. All of these results show agreement. These comparisons for the 24 inch 
test section indicate that the tangential grid scheme provides a useful grid which is 
validated by the experimental data. It should be noted that even though both of the 
grid schemes tried here satisfied the minimum desired yt value of 1 only the 
tangential grid scheme provided good agreement with experimental data. In this 
case there was little difference among grids which had yt values less than 5. 


data ysf y* 





X, In. 


(a) 201 grid points in y direction. 
Figure 74. Boundary-layer thickness as a function of longitudinal tunnel location 
computed by the ins2d for 24" tunnel at x = 16" using the tangential grid. 
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data ysf y* 





6, in. 
xX, in. 
(b) 401 grid points in y direction. 
Figure 74. Concluded. 
1/7 power law : 
— — — ins2d, y sf=32 & j=201 | 
O | : 
y/d 





(a) Tangential scale factor of 32 and 201 grid lines vertically. 
Figure 75. Boundary-layer profile computed by the ins2d for 24" tunnel at x = 16° 
using the tangential-transformation grid. 
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1/7 power law : 
— — — ins2d, y sf=32 & j = 401 


y/6 
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(b) Tangential scale factor of 32 and 401 grid lines vertically. 
Figure 75. Continued. 


1/7 power law 
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(c) Tangential scale factor of 64 and 201 grid lines vertically. 
Figure 75. Concluded. 
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b. 120 Inch Test Section 


The grid was used to represent the 120 inch wind tunnel test section 
(C8000) from the floor to the ceiling and between the longitudinal stations from 0 
to 120 inches. The grid clustering for Roberts’ transformation 1s determined by the 
value used for the stretching parameter B (equation 36). For the 120" test section, 
the boundary layer thickness (5) ranged from 1.1 to 1.5 inches depending on the 
longitudinal survey location. This leads to a recommended stretching parameter of 
about 1.05. In the ins2d code the resultant grid yields a minimum yt value of 7.0, 
which is larger than the recommended value of about 1. Using 201 grid points in 
the y direction, normal to the freestream, three values of B were used which 
resulted in a reduction in the minimum yt from 8.6 to 0.7. None of these cases 
resulted in converged calculations even after 10,000 or more iterations. As 
discussed earlier, this is an undesirable feature of these grids. The following table 
outlines Roberts’ transformation cases which were used to evaluate these grids. 


aneeneen Mae 
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The computed boundary layer thicknesses (6) are presented in Figure 
76 along with the experimental data and the 1/7 power law result. The 
experimental values are much less than those from the 1/7 power law, as noted 
earlier in the discussion of the experimental data. The computational result from 
the coarsest grid produced thicknesses which were less than one-half of that 
measured in the experiment. The computational result from the two finest grids 
where yt is near 1 provide thicknesses which are less than one-third of that 
measured in the experiment. These inadequate results reflect the lack of 
convergence in these computations and suggest that the Robert's transformation 
grids are not satisfactory for the 120 inch test section. 
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The tangential-grid scheme was also used for the 120-inch test- 
section tunnel (C8000). Scale factors from 8 to 72 were selected to achieve a 
range of yt values which approached 1. The number of y grid points ranged from 
201 to 601 to evaluate grid-density effects. In contrast to the tangential grid 
transformation for the 24 inch test section, not all of these cases converged in less 
than 100 iterations. The yt and iterations associated with the tangential-grid 
scheme are presented in the following table: 











Tangential grid C8000 y grid points 


scale factor + / iterations 


em i i 
| | savesc0ne | saviosisne | 
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pn | avse | once | 


Here nc means ‘not converged solution’ and c means ‘converged solution’. It is 
noted that while all of the converged cases required less than 100 iterations, the 
unconverged cases were run for about 8000 or more iterations. It is noted that 
with grids using scale factors of 8 and 16, the solutions did not converge 
Increasing the grids from 401 to 601 y grid points (using a scale factor of 64) 
caused the calculations to diverge even though the yt value was reduced from 0.8 
to 0.6. 







Calculated boundary-layer thicknesses are compared in Figure 77. 
Figure 77(a) presents the results for 201 y grid points. Only the unconverged scale 
factor of 8 results differ significantly from the experimental data. Even the 
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unconverged scale factor 16 results provide good agreement with the experimental 
data. As a result, good agreement is achieved with yt of 6.5 down to 1.4. Figure 
77(b) presents the results for 401 y grid points. None of the unconverged results 
agree with the experimental data. The scale factor 8 result overestimates the 
boundary layer thickness and the scale factors of 16 and 32 underestimate the 
experimental data. Only the two converged results for scale factors of 64 and 72 
achieve reasonable agreement with the experimental data. In this y grid size 
convergence was achieved for y* of 1.7 or less. The final comparison in Figure 
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Figure 76. Boundary-layer thickness computed by the ins2d for 120" tunnel using 


Roberts’ transformation grid, 201 y points. 
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| data ysf y* data ysf y™ 
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(a) 201 grid lines vertically. 
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Figure 77. Boundary-layer profile computed by the ins2d for 120" tunnel using 


the tangential-transformation grid. 


: data ysf y* data ysf y* 
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(b) 401 grid lines vertically. 
Figure 77. Continued. 
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- data ysf y* 
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(c) Tangential scale factor of 64. 
Figure 77. Concluded. 


77(c) for scale factor 64 grids show that increasing the number of y grid points 
from 201 to 401 to 601 reduces the y+ from 1.5 to 0.8 to 0.6 for the finest grid. 
Yet the worst result was from the finest grid, which is also the only unconverged 
result shown in this plot. 


In summary, the computational results from both of the test section 
lengths used in the experimental investigation, show that, first, a grid which 
provides a converged solution is needed before a reasonable comparison with 
experimental data can be expected. The value of y* is of less importance as 
shown by the Robert's grid where a yt of 0.7 gave poor results; yet a tangential 
grid with a yt of 7 gave good results. The appropriate grid clustering, size, and 
density are not always apparent before comparing results with experimental data. 
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2. Driven Cavity Computations 


A frequently used model! problem for cavity flow is the driven cavity. This 
cavity has three stationary walls and one wall which moves at a constant velocity. 
The velocity magnitude is determined by the Reynolds number of the calculation. 
Ghia et al [Ref. 67] conducted a very detailed analysis using a vorticity-stream- 
function formulation for the incompressible Navier-Stokes equations. Uniform 
grids were used in these calculations. For Reynolds numbers from 100 to 3200 the 
grid was 129 by 129; and for Reynolds numbers of 5000 or greater the grid was 
257 by 257. An example of their streamline flow pattern was presented in Figure 
10 for a Reynolds number of 10,000. In this paper, flows were computed for 
Reynolds numbers from 1000 to 10,000 and then the results were compared with 
those from several other investigators [Ref. 68, 95,and 96]. The Ghia u velocity 
results along a vertical line which passes through the center of the cavity compared 
with the other computations in Figure 78. The origins in Figure 78 are displaced - 
to separate the various curves. At low Reynolds numbers the boundary layers are 
very thick. As Reynolds number reaches and exceeds 5000, 6 asymptotes toward 
a converged magnitude. Away from the cavity walls, the velocities tend toward a 
linear variation especially for the high Reynolds number cases where the boundary 
layer thickness is a small fraction of the cavity depth. In most of these 
comparisons, the Ghia results are very consistent with those of the other three 
investigators. The most notable difference occurs when the Ghia results are 
compared with those from Ref. 95, which used the coarsest grid (50 x 50) of the 
four sets of results shown. 


Comparable calculations done using the ins2d incompressible Navier- 
Stokes code [Ref. 91 and 92] for Reynolds numbers of 1000 and 10,000 are — 
presented in Figure 79 with the results from Ghia et al [Ref. 67]. In Reference 91 
it was found that a clustered 81 by 81 grid agreed with the Ghia et al [Ref. 67] 
results. For the present investigation, an ins2d calculation was made using an 81 
by 81 grid with tangential spacing. The grid was generated using the cavity grid 
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Figure 78. Comparison of u velocity along vertical lines through the geometric 
center of the cavity [Ref. 67]. 
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generation code cgg.f as described previously. The u (Fig. 79(a)) and v (Fig. 
79(b)) velocity components on vertical or horizontal lines, respectively, are non- 
dimensionalized by the moving wall velocity. These comparisons show very good 
agreement for both velocity components. This calculation is consistent with 
previous ins2d results and confirms that consistent grids and boundary conditions 
are being used in the present investigation. 


y/D 
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(a) u velocities. 


Figure 79. Comparison between ins2d[Ref. 91 and 92] and Ghia et al [Ref. 67] 


Navier-Stokes computational results for Re of 1000 and 10,000. 
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(b) v velocities. 
Figure 79. Concluded. 


3. Open Cavity Flow Computations 


The ins2d results of the previous two sections provided consistent results 
for the two SJSU tunnel test section configuration’s (24" and 120") boundary 
layers and for a driven cavity problem. These data provided the confidence 
needed to progress to the computations of the wind tunnel with a cavity and with a 
test section diffuser. 


Initial computational efforts focused on the 24" test section with one cavity 
(D = 1"). The tunnel grid from the tunnel grid generation program stgg.f was 
used as a starting point. The single cavity was added using the cavity grid 
generation program cgga.f. This program used an 81 by 81 cavity grid with 
tangential spacing, using a scale factor of 16. Then the tunnel grid with a scale 
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factor of 32 vertically was augmented from 201 high by 495 long to satisfy the 
longitudinal spacing criteria after the cavity grid was introduced. This criteria uses 
a maximum grid spacing increase of no more than 20 percent between adjacent 
grid cells. The augmentation increased the grid to 201 high by 495 long for the 
basic 24" test section. The resultant experimental and computational boundary 
layer thicknesses are presented in Figure 80 for both the tunnel without a cavity 
and the tunnel with the 1" cavity. These results show a small increase (= 0.03") in 
6 when the cavity was present. In contrast the experimental data show a crossover 
in 6, so that at x = 5 the basic tunnel has a slightly thinner 6, while at x = 16 the 
basic tunnel has an increased 6. In the region immediately downstream of the 
cavity there is only a small difference in the 6 for these two configurations. The 
boundary layer profiles were very similar to the one shown in Figure 75(a). 


Similar computations were done for the 120" test section both with and 
without the 7° test section diffuser. In this case the same 81 by 81 tangentially- 
spaced cavity grid was used. The freestream grid was increased from 321 to 561 
by cgg.f to satisfy grid spacing change rate requirements. The ins2d code did not 
converge to a solution for this case. The freestream grid was split into two grids 
(201 by 321 each) and merged using the PEGSUS overset scheme [Ref. 94]. The 
first grid extended from x = 0" to 96" and the second grid from x = 96" to 120". 
This provided a denser grid in the test section region where the boundary layer 
thicknesses were measured and where the cavity and test section diffuser were 
located. For the cavity cases, the cavity grid was overset as a third grid using 
PEGSUS. These computations converged in about 200 to 240 iterations 
depending on the specific configuration. 


The experimental and computational boundary layer thickness results are 
presented for the test section region in Figure 81. The experimental boundary- 
layer thickness was slightly greater with the 1" cavity than without the cavity. The 
computed 6's are nearly the same as the cavity case where the experimental data 
are a small increment (= 0.09") thicker. For the case with the test section diffuser 
(Figure 82) the experimental data for the cavity case has a slightly thinner 6 (Ad 
= ().03") while the computational data has a little larger difference (Ad = 0.05"). 
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While the computational results shown in Figure 80 through 82 have some 
disagreement with a few experimental data points, the trends tend to be consistent 
with the experimental data trends for the zero pressure gradient in both test 
sections; the presence of the cavity slightly increases the boundary layer thickness. 
For the 120" test-section case with the 7° test-section diffuser, the cavity slightly 
decreases the boundary-layer thickness 6. Overall, both experimental and 
computational results show only small changes in the boundary layer 
characteristics due to the presence of a 1" cavity. 
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Figure 80. Boundary-layer thickness as a function of longitudinal tunnel location 
computed by ins2d for 24" tunnel. 
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Figure 81. Boundary-layer thickness as a function of longitudinal tunnel location 
computed by ins2d for 120” tunnel. 
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Figure 82. Boundary-layer thickness as a function of longitudinal tunnel location 
computed by ins2d for 120" tunnel with the 7° test section diffuser. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


One airfoil high-lift concept uses multiple cavities on its upper surface to 
increase camber. The usefulness of cavities in an adverse pressure gradient is 
determined by their effect on the viscous flow in the boundary layer downstream 
of the cavities. Maintenance of attached flow requires that the flow over the 
cavities and downstream of them retains enough momentum to overcome the 
kinetic energy loss due to an adverse pressure gradient, shear-layer flow gradients, 
and viscous dissipation. For this concept's feasibility to be demonstrated the effect 
of cavities on boundary-layer characteristics especially in an adverse pressure 
gradient must be determined. The present investigation concentrated on the effect 
of cavity flow on the attached surface boundary-layer characteristics in adverse 
pressure gradients. 


The purpose of the present investigation was to experimentally and 
computationally determine the effect of 1, 2, or 4 two-dimensional spanwise 
cavities with square cross-sections on the boundary-layer characteristics in both a 
constant-pressure flow and for two adverse pressure-gradient flows. The 
experimental investigation was conducted in the San Jose State University (SJSU) | 
12" by 12" low-speed tunnel. Two test-section lengths, 24" and 120", were used to 
obtain a boundary-layer thickness relative to the cavity depth 6/D which was 
either less than 0.5 ( 6/D < 0.5 ), where unsteady flow oscillations may exist in the 
cavity (or cavities), or greater than 1.0 ( 6/D > 1.0 ) where there is predominately a 
steady, standing vortex flow in the cavity (or cavities). Adverse pressure gradients 
(dp/dx > 0) were obtained using 7° and 14° inserts on the floor of the test section. 
Total pressure profiles in the boundary layer were measured at several longitudinal 
locations both upstream and downstream of the cavities. Computational results 
were obtained using a numerical solution to the incompressible Navier-Stokes 
equations. Additional tests were conducted with a single cavity which was either 
open or closed to determine the effects of surface discontinuities on the boundary- 
layer development. 


The measured profiles demonstrated boundary layer similarity for the basic 
test section and for the 7° test-section diffuser. There was a significant boundary- 
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layer profile variation from the 1/7 power law for the 14° test-section diffuser 
which was typical of an adverse pressure gradient. The single and double cavities 
showed only a small deviation from the profile measured with no cavity. The 
four-cavity configuration profile, in most cases, represented a boundary layer 
profile similar to those found in an adverse pressure gradient. All of the surface 
pressure data suggest that the presence of the cavities had no noticeable effect. It 
may be concluded that cavities may be deployed with only a small change to the 
boundary-layer profile and without significantly modifying the resultant pressure 
distribution. This important conclusion shows that this high-lift concept may be 
feasible because the multiple cavities make only small changes in the boundary- 
layer characteristics. 


Another aspect of this application is the effect of surface non-uniformities 
which could be caused by cavity doors on the boundary layer and on the surface 
pressure distribution. As expected, a poor fit of a cavity door could increase 
boundary-layer thickness and the skin-friction drag. There were differences in the 
longitudinal pressure distribution for both test sections. These results indicated 
that the surface uniformity for a closed cavity had an effect on the freestream 
characteristics, which determine the surface pressure distribution. It was also 


found that the open, cavities produced less change in the surface pressure | 


distribution than did surface non-uniformity in an otherwise plain ceiling. This is 
a further indication of the feasibility of the use of cavities. 


In summary, the computational results for both of the test-section lengths 
used in the experimental investigation, show that a grid which provides a 
-converged solution is needed before a reasonable comparison with experimental 
data can be expected. A small value of y* is of less importance than convergence. 
This was shown by the unconverged Roberts’ transformation grid with a yt of 0.9 
which gave poor results, while a converged tangential grid with a y* of 7 gave 
good results. Computations for the 120" test section with either a one inch cavity 
and/or with a 7° test-section diffuser showed only small changes in the boundary- 
layer characteristics. The appropriate grid clustering, size, and density were not 
always apparent before comparing computational results with experimental! data. 
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Stated another way, it was found that the experimental data were needed to guide 
the computational effort. 


This experimental and computational investigation has shown that multiple 
open cavities in either zero or adverse pressure gradients make only small changes 
in the boundary-layer characteristics. It was also found that a small backward- 
facing step made larger, adverse changes to the boundary layer than those made by 
the cavities. It is concluded that multiple cavities may be a feasible high-lift 
concept. It is recommended that the cavity concept be implemented in an airfoil 
and that it be tested in a wind tunnel to quantify its high-lift and drag 
characteristics. 
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APPENDIX A 


TEST RUN SCHEDULE 


The run schedule for the experimental investigation is presented in this 
appendix. It lists the following information: 


Run Sequential number assigned to a group of data 
Configuration A four digit code which describes the tunnel configuration 
code The first digit identifies the test section length: 


0 or | refers to the 24" long test section 
8 refers to the 120" long test section 
The second digit identifies the diffuser angle: 
O means basic test-section with no diffuser 
7 means the 7° test section diffuser 
14 means the 14° test section diffuser 
The third digit identifies the cavity depth in inches: 
OQ means no cavity 
1 means 1 inch deep cavity(-ies) 
2 means 2 inch deep cavity 
3 means 2 inch deep cavity with a movable door which 
can be used to close the cavity or open the cavity 
to the freestream flow (i.e., c-closed or 0-open). 
The fourth digit identifies the number of cavities in the 
streamwise direction as either 0, 1, 2, or 4 cavities. 


199 








Run type This identifies the primary data taken as either CP, surface 
Cp distribution, or BL, boundary layer profiles. 


boundary layer This identifies the x location of the boundary layer probe tip. 


probe position 

b.l. probe height = This identifies the vertical dimension of the boundary layer 
probe. | | 

b.l. rake installed This identifies whether the b.1. rake is installed at the x=16" 
location. 

Notes Notes related to the run describing unexpected items or 


information unique to a particular run which is not described 
in the other columns. 


An entry is made in the run schedule only when there is a change in a 
parameter. If there is no entry then the parameter is unchanged from the previous 
run. 
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code Type layer probe head installed | | 7 
nosition height | 

ee Rs 
PO of OOOO. CR INI IN i ee! 
fT 102 T s—C‘iE:C“C( CB void run 102, invalid inputcode 
OS eee eee as ee a 
Pe ON sc et ae ee ee 
FO ee rll na ac I ge Mc 
OO ees fe IO Nc ee a ee 
Pos [0700 sf CP N/A NVA [i void run, scanivalve did not step properly | 
C109 fT t—(tiséitéisCSi‘(‘“‘(C:é‘COCTO TT repeated run with correct scanivalve stepping | 
Po tC BE 6 0.022 [si pressure sample time increased to 20 sec.to_ 
Pill fC ef  <cont'd> improve scanivalve data repeatibility | 
ew seca fe Ne 
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Bad cavity data-used 20 sec. sample time-void run 


Scanivalve sample time reduced from 8sec to 2sec 


Repeated with tape removed from ports 1 & 16 


epeated with transducers Pt & Pt,bl exchanged 


R 
Pt,bl tube replaced before run. Leak detected. 
Ps,bl tube replaced before run. Moisture detected. 
Repeated run 151 w/ new Ps,bl tube. 

Repeated run 149 w/ new Ps.,bl tube. 

Repeated run 148 w/ new Ps,bl tube. 


Repeated run 147 


sed wrong configuration code - 0014 
epeated run 159 
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Run 


Configuration 
code 


Run |} 
Type 


boundary 
layer probe 
position 


b.1. probe 
~ head 
height 


b.]. rake 
installed 





Notes 





Pt,bl leak. 


epeated run 159; lost some Cp data. 
epeated run 156; lost data. 

epeat run 186 

epeat run 186 


re) Ve) 


epeatibility check of run 189. 
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Repeat run 215. 
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epeat of run 222 w/ BL 


probe @ x= 


108 in. 
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Run Configuration Run boundary 
| code Type layer probe 
position 


“8014 
ee eee eens ee 
230 fC BE 


b.l. probe 
head 
height 


b.L. rake’ 
installed 


Notes 


qeo started to vary +0.2 psfg. Void data. 
epeat run 222 w/ BL probe @ x=108". Tape 
cont'd>obstruction ahead of the probe in run 228 


A| A 


epeat run 234 


his run used to calibrate new Pt,bl probe 


Repeat run 240 w/ new Pt,bl probe. 
Repeat run 229 w/ new Pt,bl probe. 


Void run, used the wrong config. code-8014. 
Repeat run 247 with right config. code-8000 


NOTE: Void runs 249-257 due to bad Pt,bl probe. 


Repeat run 252 with new Pt,bl probe. 

Repeat run 253 

New Pt,bl probe w/ 0.034" probe height used. 
Tried a modified 0.022" probe 

Tried a second 0.022" probe 
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layer probe head 
Hosition | height 


alid data w/ new Pt,bl probe - 0.034" height 
oid data, probe deflected against the ceiling. 
epeat run 263, probe just touching the ceiling. 
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oid run, delta > 2" 

epeated run 285 

epeated run 285 

epeated run 284, checks out o.k. 
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29s | TCTC~wmSC“‘(‘SSCrLCOOCOC*d:CNW*C‘CW*‘LCU CNC Repaired bl probe and obtained valid data__| 
306 8a Pda dT _ | Void rn because of wrong configuration code 
997 so” [Repeated 296 
door closed: fe fg 
dooropen [| CP {| 
RE [02 ee Caan Ee =) a (EAE UE: SNE! | Nee eee, 
Sn (1) CR (meena) (SNe (TERNS) (Ne AA UN (Een: 
[302__| doorclosed [|S CC~dEC‘(‘(SCL OC Odd reflex in bl survey for O.l’<y<0.37 

5 (203 Adjusted the trap door to flush & repeated run 302 
S | 304 [ 843ic cp ia Repeatedrun 297 

[305 | C™C™~—“—CCSTSCéSCdTs TC~sdEC $ST Repeatedrun 298 
Oller ee ON a 
Sloper ee ee 
Fee cee OO i al a ce 
[309 [| doorclosed [OT CdS (ST 
T 310 | | CP [112 [  [ {Repeated run 297 after repair of tube @ x=11.03" 
SiS __| Repeated run 310 w/ I probe flush to ceiling 
[312 [| dooropen [| SCC Cd Cd Repecatedrun 299 
oi ee ee 
i 
3s ff té“<~*sSC‘“‘CSC*LSC*C‘LOBC Tb profile had odd reflex as seen in run 302 
T 3i6 of t~“‘;SCS*rTSC‘“‘(NSNNC*SE‘C(N$’$#$NS~=;WUNCSC‘#$§$WM[/ ws sf Repaired loose tape & repeated run 297 
Co 
38 os OO 
319 | TP C—wSCOCSsSSCSCSCSCSCSTTT....C ''"'"'|__ 
0 ee 
eee ee 








80C 


boundary | b.l. probe 


Run Configuration 
code layer probe head 
position height 


Gol Gof Gol Gol Gol Go} Go 
Cn} Cat Cn} By Be BY 
bol} CE NOT OO] ~ I ON 


se) 
Nn 
ee) 


aE! 10075 [ 
SE 1S AS 
CE NGS CNRS SARIN 
a SP 
09235 
Se A 7 
I VO ae SAAR 
05 
18 
ff 0 
109225 
[80sec 
[door open_| | id 
ee Ge: a TE: RE 
SS 
dooropen [| Cd 
SE 
doorclosed [—_ [dP 
(DS 
[door open || sd 
[_o03io [| _cP_| 12 
door closed_|BL_|____| 
i Ma Te SANE 
[door open [i Sd 
a NS (: 
doorclosed || sd 
a A 
dooropen_ [| sd 
iecccacaeSe SORS HENGE:) *- sel RA 
(door closed [ff 





b.1. rake 
installed 
no 





P 


itot static tube @ y= 


ore 





60C 





Run Configuration Run boundary | b.l. probe] b.l. rake Notes 
code Type layer probe head | installed 
position height 


354 0000 16 tunnel diffuser w/o bl rake fairing 





2 >) ae Traversed bl probe to y=9"; pitot @ y=4.5" 
Ee cee ee 
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nn ae 5.5 Traversed to y=9" 
359 16 
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Used coarse traverse to y=9" 
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Use extended y traverse 
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Extended y traverse. 
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APPENDIX B 


SURFACE PRESSURE TAP LOCATIONS 


The ceiling surface and cavity pressure tap locations used to obtain pressure 
distributions in the present investigation are listed in this appendix. It includes the 
following tables: | 


















number | number Surface 
of 


ceiling 






of cavity 






taps 






N/A 
| 16 | 8 | one 2" deep cavity w/ trap door and adjacent surfaces | 
, 16 44 


Pe | four 1" deep cavities and adjacent surfaces 
The basic test section of this wind tunnel uses only the basic ceiling whose pressure 
tap locations are given in Table 3 and a plain uninstrumented tunnel floor. When 
an adverse pressure gradient is desired one of the test section diffusers whose 
pressure tap locations are given in either Table 1 or Table 2 is added on the tunnel 
floor along with a modified contraction at the test section entrance. When a cavity 
configuration is desired the basic ceiling is replaced by the appropriate ceiling 
given in Tables 4 through 8. All of the x station dimensions in these Tables are 
given for the 24" long test section. The 120" long test section is achieved by 
adding a 96" insert with a variable floor location upstream of the original test 
section. As a consequence the x stations for the 120" long test section may be 
obtained by adding 96" to the x station values given in these Tables. 
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Table 1. 7° test section diffuser pressure tap locations: 


x, inches y; inches 





16 


Table 2. 14° test section diffuser pressure tap locations: 


Tap 


_———— — 
8982 
3. or. $< 
rs 
























16 21 504 -5. 381 
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Table 3. Basic test section ceiling. 


: x, Inches , inches 
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Table 5. One 2" deep cavity with a trapdoor and the adjacent ceiling surface. 


Sons Tap | x, inches , Inches | Cavi ‘Tap. x, inches , inches 
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Table 7. Two 1" deep cavities and the adjacent ceiling surface. 








“Tap. 
Ta Tap number 
| it | 1969 | 6 | ft | 6 {| 647 | 1 
2 | 2969 | 6 2 BST 
| 3 | 3.9388 | 6 8 T6938 
tE 4 | 4969 | 6 4 622 693 
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Table 8. Four 1" deep cavities and the adjacent ceiling surface. 
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APPENDIX C 


SURFACE PRESSURES 


The pressure distributions measured on the wind tunnel surfaces in the 
cavities and on the adjacent freestream surfaces are presented in this appendix. 
The data from each cavity are presented in three contiguous plots: (1) front cavity 
wall; (2) cavity floor; and (3) rear cavity wall. The freestream surface plots are 
presented for all of the pressures measured on the y = 6 inches surface. The data 
are faired only for the surfaces upstream and downsteam of the cavity (or cavities). 
The following table relates the configurations and their descriptions with the figure 
numbers: 


Test Number of Surface Cavity 
Section Cavities Pressures | Pressures 
Length | 





pale 













Test Number of 


Section 


Config- 
uration 


Length | | 
lin, | C21 C22 
1in. C32 


8431 2 in. C49 C50 


Surface Cavity 








Cavities 





Pressures Pressures 
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Figure C1. 





_—t— Upstream 





Pressure distribution on the tunnel surfaces adjacent to the 0011 
cavity. 


Cavity 


-0.08 





4 (6.938 in.) 7 6 
(6.50in.) 6.5 7.5 (7.50 in.) 


MKS 


Figure C2. Pressure distribution in the 0011 cavity. 
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Figure C3. Pressure distribution on the tunnel surfaces adjacent to the 0711 
cavity. 


Cavity 





Y: 6 7 (6.938 in.) 7 6 
x (6.50in.) 6.5 7.5 (7.50 in.) 


Figure C4. Pressure distribution in the 0711 cavity. 
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Figure C5. Pressure distribution on the tunnel surfaces adjacent to the 1411 
cavity. 
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Y: 6 7 (6.938 in.) 7 6 
x: (6.50in.) 6.5 7.5 (7.50 in.) 


Figure C6. Pressure distribution in the 1411 cavity. © 
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Figure C7. Pressure distribution on the tunnel surfaces adjacent to the 8011 
cavity. 
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Figure C8. Pressure distribution in the 8011 cavity. 
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Figure C9. Pressure distribution on the tunnel surfaces adjacent to the 8711 
cavity. 
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Y: 6 7 (6.938 in.) 7 6 
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Figure C10. Pressure distribution in the 8711 cavity. 
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Figure C11. Pressure distribution on the tunnel surfaces adjacent to the 8411 


cavity. 


Cavity 





Y: 6 7 (6.938 in.) 7 | 6 
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Figure C12. Pressure distribution in the 8411 cavity. 
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Figure C13. Pressure distribution on the tunnel surfaces adjacent to the 0012 
cavities. 
Cavity #1 Cavity #2 
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Figure C14. Pressure distribution in the 0012 cavities. 
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Figure C15. Pressure distribution on the tunnel surfaces adjacent to the 0014 
_ Cavities. | 
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(a) Pressure distribution in the first two cavities. 
Figure C16. Pressure distribution in the 0014 cavities. 
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(b) Pressure distribution in the second two cavities. 
. Figure C16. Concluded. 
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Figure C17. Pressure distribution on the tunnel surfaces adjacent to the 0714 
cavities. 
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Cavity #1 Cavity #2 





(a) Pressure distribution in the first two cavities. 
Figure C18. Pressure distribution in the 0714 cavities. 


Cavity #3 Cavity #4 





(b) Pressure distribution in the second two cavities. 
Figure C18. Concluded. 
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Figure C19. Pressure distribution on the tunnel surfaces adjacent to the 1412 
cavities. 
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Figure C20. Pressure distribution in the 1412 cavities. 
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Figure C21. Pressure distribution on the tunnel surfaces adjacent to the 1414 — 
cavities. | | 


Cavity #2 





(a) Pressure distribution in the first two cavities. 
Figure C22. Pressure distribution in the 1414 cavities. 
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Cavity #3 Cavity #4 





(b) Pressure distribution in the second two cavities. 
Figure C22. Concluded. 
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Figure C23. Pressure distribution on the tunnel surfaces adjacent to the 8012 
cavities. 
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Figure C24. Pressure distribution in the 8012 cavities. 
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Figure C25. Pressure distribution on the tunnel surfaces adjacent to the 8014 
cavities. 
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(a) Pressure distribution in the first two cavities. 
Figure C26. Pressure distribution in the 8014 cavities. 
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(b) Pressure distribution in the second two cavities. 
Figure C26. Concluded. 
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Figure C27. Pressure distribution on the tunnel surfaces adjacent to the 8714 
: cavities. 
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(a) Pressure distribution in the first two cavities. 
Figure C28. Pressure distribution in the 8714 cavities. 
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(b) Pressure distribution in the second two cavities. 
Figure C28. Concluded. 
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Figure C29. Pressure distribution on the tunnel surfaces adjacent to the 8412 
cavities. 
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Figure C30. Pressure distribution in the 8412 cavities. 
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Figure C31. Pressure distribution on the tunnel surfaces adjacent to the 8414 
cavities. 
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(a) Pressure distribution in the first two cavities. 
Figure C32. Pressure distribution in the 8414 cavities. 
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‘(b) Pressure distribution in the second two cavities. 
Figure C32. Concluded. 
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Figure C33. Pressure distribution on the tunnel surfaces adjacent to the 0021 
cavity. 
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Figure C34. Pressure distribution in the 0021 cavity. 


238 








Figure C35. Pressure distribution on the tunnel surfaces adjacent to the 8021 
cavity. 


Cavity 





Y: 6 g (8.00 in.) g 6 
X: (103.5 in.) 103.5 105.5 (9.437 in.) 


Figure C36. Pressure distribution in the 8021 cavity. 
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Figure C37. Pressure distribution on the tunnel surfaces adjacent to the 8421 
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Figure C38. Pressure distribution in the 8421 cavity. 
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Figure C39. Pressure distribution on the tunnel surfaces adjacent to the 1421 
cavity. 
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Figure C40. Pressure distribution in the 1421 cavity. 
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Figure C41. Pressure distribution on the tunnel surfaces adjacent to the 0031 
cavity. 
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Figure C42. Pressure distribution in the 0031 cavity. 
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Figure C43. Pressure distribution on the tunnel surfaces adjacent to the 0731 
cavity. 
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Figure C44. Pressure distribution in the 00731 cavity. 


243 















0.5 , , , 
——O—- Upstream, open : 
——}-—— Downstream, open: 

0.4 6 = Upstream, closed ae Ditcees "eile ouenaaemacceoaal 
— @ — Downstream, closed i 

0.3 ena en edn MERC” cARNU RNIN RSE Rint rs cr meee ee 

C 
p 
0.2 
0.1 
0 
0 5 10 15 20 25 
X, in. 


Figure C45. Pressure distribution on the tunnel surfaces adjacent to the 1431 
cavity. 
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Figure C46. Pressure distribution in the 1431 cavity. 
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Figure C47. Pressure distribution on the tunnel surfaces adjacent to the 8031 
cavity. 
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Figure C48. Pressure distribution in the 8031 cavity. 
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Figure C49. Pressure distribution on the tunnel surfaces adjacent to the 8431 
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Figure C50. Pressure distribution in the 8431 cavity. 
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APPENDIX D 


BOUNDARY LAYER PROFILES 


The boundary layer data from this investigation is presented for each wind 
tunnel configuration in this appendix. There were four to six longitudinal 
locations for each configuration where boundary layer profiles were measured. 
Each figure identifies the x location for the data presented. The displacement 
thickness was computed from the data and then used to non-dimensionlize the 
distance y which was measured from the wind tunnel ceiling. A pitot probe 
measured the total-pressure variation and the static pressure was obtained from a 
port located on the ceiling at the appropriate longitudinal location. As expected, 
the static pressures were nearly constant and the total pressures varied from a few 
pounds per square foot greater than the static pressure near the surface to the 
freestream total pressure at the boundary layer edge. The following table relates 
the configurations and their description with the figure numbers: 


Config- Test Number of 
uration Section Cavities 
Length 
8000 


8400 120 in. 
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OO11 


0021 


0711 


1411 


1421 


8011 


8021 


8711 


B41 


8421 


0012 


0014 


0714 


1412 


1414 


8012 


8014 


s) 
ct 
S 
5 


8714 


Cavities | e 
el“ 
ef te 


ho 
aa 
p_i © 
RP 
ttt 


, 

Cee toe Te 
vin | 1 | sm | 
om | 1 | am |p 
in | 2 | in | 
min | | im | 
pin | as | in | 
min | | im | 
min | a | im | 
| voin | 2 | sm | » 
von | 4 | in | 
von | | im | 
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Notes 

























No. of 
Cavities 






Config- 
uration 
Se 
Lisa | on | a | 
coste | min | 1 | 2 
[Loose | ain | a | 2 
Lome | wm | 1 | 2 
oa |: 
a 
1 |. 
os 


Cavity Figure No. Notes 
Depth 













07310 24 in. 
1431c 24 in. 
1431o 24 in. 
| 8031¢c 120 in. 
80310 120 in. 


8431c 120 in. 
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Figure D1. Boundary layer profiles for 0000 configuration. 
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Figure D4. Boundary layer profiles for 8700 configuration. 
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Figure D6. Boundary layer profiles for 8400 configuration. 
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Figure D8. Boundary layer profiles for 0021 configuration. | 
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Figure D10. Boundary layer profiles for 1411 configuration. 
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Figure D11. Boundary layer profiles for 1421 configuration. 
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Figure D12. Boundary layer profiles for 8011 configuration. 
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Figure D14. Boundary layer profiles for 8711 configuration. 
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Figure D16. Boundary layer profiles for 8421 configuration. 
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Figure D17. Boundary layer profiles for 0012 configuration. 
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Figure D18. Boundary layer profiles for 0014 configuration. 
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Figure D20. Boundary layer profiles for 1412 configuration. 
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Figure D22. Boundary layer profiles for 8012 configuration. 
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Figure D23. Boundary layer profiles for 8014 configuration. 
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Figure D24. Boundary layer profiles for 8714 configuration. 
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Figure D26. Boundary layer profiles for 8414 configuration. 
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Figure D28. Boundary layer profiles for 00310 configuration. 
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Figure D30. Boundary layer profiles for 07310 configuration. 
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Figure D32. Boundary layer profiles for 1431 o0configuration. 
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Figure D34. Boundary layer profiles for 80310 configuration. 
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Figure D36. Boundary layer profiles for 84310 configuration. 
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